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The theory of the interpretation of seismic travel-time curves for refracted rays 
in horizontal structures is treated after the manner of Herglotz-Wiechert, under the 
customary assumption that the ray paths obey the laws of geometrical optics. Mul- 
tiple valued travel-time curves, discontinuous velocity functions, and the discon- 
tinuous travel-time curves associated with a slower speed bed receive special consider- 
ation. It appears that interpretations satisfactory from the theoretical point of view 
may be obtained in these cases, although, experimentally, sufficiently complete data 
to meet the requirements of theory may often be difficult or impossible to obtain. 


I* 1907, the problem of the direct interpretation of seismic travel-time 
curves for a symmetrical earth was treated by Herglotz,! who perceived 
that the determination of the velocity of elastic waves in the earth's interior 


1G. Herglotz, Phys. Zeits. 8, 145 (1907). 
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could be identified with the solution of Abel’s integral equation. In 1910, 
Bateman’ independently solved the same problem, and by essentially the 
same method. Finally, in 1910, Wiechert* brought the subject essentially to 
its present state, by a contribution of a practical nature. He succeeded in re- 
ducing the integrations involved which had been very tedious, to a quite 
simple and elegant process. 

The special nature of the problem which Herglotz wished to solve per- 
mitted him to restrict his analysis in certain respects; first, the velocity was 
supposed always to increase with increasing depth; second, the variation of 
velocity with depth was assumed to be of such a type that no ray paths inter- 
sected. While these assumptions were justifiable in Herglotz’s problem, there 
are other cases of the same general problem for which such restrictions may 
not be imposed. For example, in seismic prospecting, the presence of a slow 
speed layer leads to violation of the first requirement, and an acceleration in 
the velocity-depth function may lead to the violation of the second. Again, in 
the recent explorations of the upper atmosphere by acoustic travel-time data, 
and in the study of layering in the upper fifty miles of the earth’s crust by 
means of local earthquakes, the two conditions mentioned need not in general 
be satisfied. In closing his 1907 paper, Herglotz expressed his intention of 
examining the problem further, free from these restrictions. However, a 
search of the literature seems to show that little further attention has been 
given to the theory of the subject. 

In discussing the two restrictions noted above, it will be convenient to re- 
cast the problem slightly, by assuming a flat instead of a spherical earth. For 
the newer applications, such as those in seismic prospecting, in the seismic 
studies of the upper crust, and in the acoustic study of the upper atmosphere, 
it is usually appropriate to neglect the curvature of the earth. Corresponding 
to Herglotz’s assumption of spherical symmetry, the assumption of horizontal 
uniformity is made, so that in the present treatment the velocity is permitted 
to vary only with depth. 

It will be convenient to outline briefly the well-known method of Herglotz- 
Bateman-Wiechert, in order that some of the relationships that occur may 
be conveniently available for subsequent reference. 


Part I. NORMAL TyPpE oF TIME-DISTANCE CURVE 
1.1 Review of Herglotz-Wiechert solution 


The problem treated by Herglotz may be visualized by reference to Fig. 1. 
The source or focus is at the ground’s surface. The velocity of the elastic 
wave is supposed to increase with depth, as indicated by the graph in the 
lower left. As in the optical analogue, the energy is assumed to be propagated 
along the curved ray paths shown, in accordance with Fermat'’s principle ot 
least time; and for the present case of a horizontally uniform medium, it is 


? H. Bateman, Phil. Mag. 19, 576 (1910). 
5 Wiechert and Geiger, Phys. Zeits. 11, 294 (1910). 
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easily shown that Fermat’s principle implies Snell’s familiar law of optical 
refraction; that is, along any given ray the condition 


sin i/v = sin ip/% (1) 


is satisfied, where v is the velocity at the surface, and i) the initial angle of 
incidence, or, by symmetry, the angle of incidence at the point of emergence 
of the ray. The value of sin 7 as a function of epicentral distance is easily ob- 
tained from the travel-time curve. For, in Fig. 2, if S represents the surface 
of the wave front just before emergence, the ratio of the wave velocity 
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Fig. 1. Normal type of travel-time Fig. 2. 


curve. Normal path. 


vp to the apparent velocity of propagation along the surface, 3, as read from 
the travel-time curve, is sin Zp. 
09/0 = SiN io. (2) 
On substituting (2) in (1) we obtain the following relation characterizing any 
given ray: 
sin i = v/d. (3) 
Thus at the deepest point of any path, where sin i=1, the wave velocity 
is precisely the travel-time velocity at the point of emergence of this ray. It 
remains, therefore, only to determine the maximum depth, z,, attained by a 
ray of travel-time velocity 3,. 
The horizontal distance, x, traversed along any ray from the origin to a 


depth gz, is 
a | 
t= f tan idz 
0 


which, by using Eq. (3), may be written 
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- f- ° (4) 
TJS Bli/et — 1/0} 


When the upper limit z, is taken as the deepest point of the ray, this equation 
may be identified with Abel’s integral equation, 


f(y) = f “ulr)ar/(y —r)".,0<n<l (A) 


whose solution is* 


sin ur d é y)dy 
u(t) = ee _Sy)dy 


r dtJe (t—y)"_ (B) 


If we write A= 2vx,, so that A is the emergence distance of the ray of apparent 
velocity 5, and introduce the substitutions 


p=1/v?; § = 1/0; f(g) = (4/2)-9, (S) 


then Eq. (4) becomes 
&  (dz/dp)dp 
fe) = f a (6) 
I/vg (p — £) 


where p is a function of z alone, and £ a function only of A. Abel’s solution of 
(6) is: 





5 1dz 1 =f" f(&)dé 
— t dp m® dpJiye (op — §)' 
or 


7 
1? — fede " 


— ® Saeed (E — p)"? 


the constant of integration being zero. The quantities which enter the integral 
in Eq. (7) are all known, and hence this equation is the desired relation ex- 
pressing z as a function of p, and hence also of v. However, the numerical 
evaluation of the integral was tedious until Wiechert*® pointed out the follow- 
ing convenient method. 

Eq. (7) is rewritten in terms of the original notation of Eq. (5) above and 
then the velocity, v,, corresponding to a given value of the upper limit p is 
replaced throughout by the equal apparent velocity, 0,. Then, noting that 
bd(1/6)* =2d(1/d), Eq. (7) becomes 


a 





1 1/vp Ad(1/d) 
Pe Jaye, [(1/8)* — (1/6,)2}* 
or, 
1 20 
Zp=- — | , A d[cosh— 5,/i]|. (8) 
T /# cosh "Ep/%, 


* Maxime Béocher, An Introduction to the Study of Integral Equations, Cambridge Univ. 
Press, 1909, p. 9. 
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If this is integrated by parts the integrated term vanishes, since A=0 when 
}=vo. Hence 


Ap 
Zp = (1/n) f cosh~! (5,/0)dA. (9) 
0 


Here A, is the epicentral distance for which the apparent velocity is i,. This 
is essentially Wiechert’s well-known formula, except that his term ro log 
[ro/(ro—z) | which occurs on the left-hand side in the case of a spherical 
earth’ has been replaced by its limit z as 79 approaches infinity.5 

The integration indicated is very simply carried out graphically, after 
reading the required inverse hyperbolic cosines from a table. Fortunately, 
relatively few points usually suffice to define the cosh (#,/v) curve within 
the needed accuracy, so that in practice the formula (9) furnishes a simple 
and rapid method of interpreting the travel-time curve. 


1.2 Validity of the solution above 


Bocher* gives as the necessary and sufficient conditions that Abel’s Eq. 
(A) above, have the continuous solution (B) in the interval a<£<b), the fol- 
lowing three: (a) f(£) must be continuous in the interval; (b) f(a) =0; (c) 
Se f(€)dé/(x—£)-" must have a continuous derivative in the interval. The 
condition (c) is fulfilled whenever (a) and (b) are true and f(£) has a deriva- 
tive which is finite with at most a finite number of discontinuities.* For 
brevity, we shall refer to these as Bécher’s conditions.® 


5 See also Maurice Ewing and L. Don Leet, Geophysical Prospecting A.I.M.E., 1932, p. 269. 

6 These conditions are, however, somewhat more severe than need be, because only con- 

tinuous functions f(£) and finite’ functions u(p) are contemplated. The function f(¢) may, how- 

ever, possess finite discontinuities without invalidating the solution, as may be easily verified 

by the following example, suggested to me by Professor O. D. Kellogg. Suppose that f(£) is 
discontinuous, as follows: 

f(}) =0 when 0O 

f(é) =1 when a 


Then, by assuming that the solution (B) holds: 


E<a 
é. 


IIA IIA 


u(p) = 0, when OSpSa 


i @ ep = 6dp 
and u(p) = —— | ——— when a Sop. 
m dp“ a(p — &)'? 


If the above values for u be substituted in Eq. (A) we get the check: 
f(t) =0 when OSE<a 


IIA 





1 p dp 
ana sce) = — [ =1 when é>a. 
mM a(p — a)(— — p)'? 


Clearly, any prescribed function f(¢) having a definite number of finite discontinuities may be 
expressed as a continuous function plus discontinuous functions of the step type used above. 
Hence the solution for the Abelian integral equation formed for such a function f(€) will be the 
sum of component solutions each obtained by the usual formula (B). 

The solution of Abel’s equation when the condition f(a) =0 is not satisfied is given by Gour- 
sat, Acta Mathematica XXVII, 131 (1903). 

7 M. Bocher, reference 4, theorem 2, p. 4. 
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Part II. A TRANSITIONAL CASE 


2.1 A focus point as a discrete point 


A second type of travel-time curve will receive mention despite its very 
special nature, because it is the transitional case between the first type men- 
tioned, and the third type presently to be discussed. In it, all the rays are 
brought to a common focus, both as regards time and distance, so that the 
time distance curve degenerates into a single point, marked P on Fig. 3a. The 
velocity depth function which produces this result is easily found on writing. 


Eq. (4) in the form 
%™m 9(dz/dv)dv 
= [oon (10) 


1, (Om? — v2)1/2 
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Fig. 3a. Travel-time curve a single point. Fig. 3b. Focus point in a continuous travel- 
Special case. time curve. 


setting x» =7o/2, a constant, v?=£ and v,,2=A; and solving the resulting in- 
tegral equation, 


f (dz/dv)dé (11) 


2-8! 


for dz/dv. It results that 


v = % cosh (12/19) (12) 


is the required function which brings all rays to the common emergence dis- 
tance A=/7p. 


The elapsed time, ¢/2, corresponding to the half range x, in Eq. (10) is 


tm dz ™ Ym-(dz/dv)dt 
ne -f oulds/de)de (13) 
Zz 0 v 


v COS 7 








4, 2(Dm? — v2)1/2 
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where z, denotes the depth to the bottom of the ray path. By using 
dz/dv=(ro/m) 1/(v? —v9?)'? from Eq. (12), and putting 7=(v/vm)*, Eq. (13) 
becomes : 





t ro 1 dn 
2 | 14. \8/8, . fa/a WALD (14) 
2 2mm SF (vo/ om) ML = n)t!*(9 — (00/%m)*)*? 


This yields, on integration, 


l To Um ro 
—< 2 -——— wT = - 
2 270m Vo 2% 
Thus the times of emergence, ¢, are also equal for all rays, and the travel-time 


curve reduces to the single point, A=ro, 4=7ro/vo. The locus of such singular 
points is, then, the straight line ¢=A/vp. 


2.2 A focus point as a point in a continuous travel-time curve 


It is of interest to examine whether a focus* point such as this can be a 
point in a continuous travel-time curve through it. In Eq. (11) the half range 
Xm Was expressed as a function of A, say W(A), where it will be recalled \=v,,2. 
(a) If v(z) be assumed continuous and increasing, dz/dv in Eq. (11) will be 
finite, and W(A) will be continuous. (b) Furthermore, dy/dd will be continu- 
ous whenever dz/dv is a continuous function of — and has a finite derivative 
which possesses at most a finite number of discontinuities.® (c) By Eq. (11), 
YW (v0") =0. 

With these assumptions, then, Bécher’s three conditions on Y(A) are satis- 
fied, so that the Eq. (11), i.e., 


r (dz/dv) dé , 
YA) = a, (\ — pa (11a) 


possesses a unique continuous solution. But none of the conditions (a), (b), 
or (c) above is invalidated by allowing ¥(A) =x,, to assume arbitrary values 
involving continuous positive or negative values of its derivative, dy/dd. In 
particular, we may assign to Y(A) a constant value over a portion of the range 
in X. But this condition on (A) is the condition that there shall be a focus 
point. We therefore conclude that Y(A) may be arbitrarily assigned subject 
to the conditions (a), (b), and (c) above to insure the existence of focus points 
at any desired epicentral distances. The requisite velocity function for estab- 
lishing this result is yielded by solution of Abel’s equation for dz/dv. 

A travel-time curve may therefore consist of two segments meeting at a 
focus point. Since v(z) was assumed increasing, the slope of the travel-time 
curve will be decreasing, and there will appear to be a discontinuity in slope 
at the focus point, as in Fig. 3b. However, in interpreting @ travel-time curve 
of this type, it is clear that the slope must be regarded as changing continu- 

§ It will be sufficient here to consider only the focusing of the rays with respect to distance, 


because as is pointed out in paragraph 3.3, a focus in time is always implied by one in distance. 
® M. Bocher, reference 4, theorem 3, p. 5. 
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ously at the focus point through all values between the two slopes on either 
side of this point. 


Part III. REVERSED TRAVEL-TIME CURVES 


Whenever the half range x (in Eq. (10)) of the rays decreases locally with 
z, crossing of the paths will occur, as illustrated in Fig. 4. This condition pro- 
duces a travel-time curve which reverses on itself as at point 1 (Fig. 4) and 
continues, always concave upward, from point 1 to point 2 where it may again 
reverse and assume the normal concave-down type. The segments meeting 
at the cusps at point 1 and 2 are always mutually tangent. The segment 1, 2 
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Fig. 4. Triple valued travel-time curve. Fig. 5. Examples of velocity function 


which produce the triple valued curve. 


will be referred to as a “reversed segment,” because rays proceeding from the 
source at steepening angles of dip describe this segment in the direction of 
decreasing A. This type of travel-time curve has been long recognized in pure 
seismology,'® but its interpretation is usually (erroneously) held to be beyond 
the scope of the Herglotz-Wiechert theory."' Its occurrence in seismic pros- 
pecting appears to be a much more common phenomenon than is generally 
supposed. 


3.1 Condition for the occurrence of a reversed segment 


The half range x traversed by a ray in attaining its vertex was expressed 
by Eq. (10). At the vertex, in order that the ray may turn back towards the 


10 B, Gutenberg, Handb. d. Geophysik, Borntraeger IV-I, 62 (1929). 
1 See reference 14. 
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surface, the velocity must obviously be increasing with depth. Hence dx/dz 
will have the same sign as dx/dv,,; therefore, the criterion for crossing may 
be taken to be dx/dv,,<0; and for non-crossing, dx/dv,>0. This derivative 


is found to have the following value on differentiating Eq. (10) and then in- 
tegrating by parts 


dx (dz/dv) °m 2¢/dv? \. 
Jon ‘ : + f } _— 


Faas = v,J. 15 
dim (Ym? ne Vo?) 1/2 -" (Vm? = y?) 1/3 ( ) 








Here (dz/dv)) denotes the slope of the depth-velocity function at the surface. 
The condition for the existence of a reversed segment is, then, that the quan- 
tity J above be negative.” 

(1) When the slope at the surface, (dz/dv)o, is negative, the ray paths will 
be concave downwards initially, and the shallowest ray will penetrate to such 
a depth that the velocity at its apex is equal to the surface velocity, vm =vp. 
With this value for v,,, the first term in the expression for J is negatively in- 
finite. Furthermore, if the velocity function has a curvature of constant 


negative sign, the integral 
*,  dz/dy? 
dv 
vy (Om? — v2)1/2 





° 


will also be negative, and the existence of a reversed segment is assured. 

(2) For positive values of the initial slope such as generally occur in the 
crust, the reverse segment can occur only if d*z/dv? is negative over part of 
the range. The values which this second derivative takes near the upper limit 
of the integral will clearly be especially influential because of the factor 
(vm? —v*)/? occurring in the denominator. 


3.2 The likelihood of the occurrence of the reversed segment in practical 
cases 


Furthermore, it should be noted that the occurrence of a layer throughout 
which the velocity is nearly constant, (d*z/dv* being assumed negative) is very 
apt to lead to crossing of the rays. For dz/dv will then be very large, and if 
R denotes the radius of curvature, 


d*z = [1 + (dz/dv)?]*/2 (2) P 
~ Ado jf 

nearly, so the absolute value of d*z/dv* will be very large under these circum- 
stances. To illustrate quantitatively the strong tendency for the occurrence 
of the triple-valued type of curve the following numerical example will serve. 

Example. Let the surface velocity serve as the unit for velocity, v=1. 
During the depth interval 0<2<1, let dz/dv be a constant, so, so that the 
velocity at unit depth is v,; = 1+1/so. Thereafter, let d*z/dv? be some negative 





dv? R 


12 This condition on the velocity function corresponds to Herglotz’s condition on the 
travel-time curve,! d*t/dA?>0. 
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constant, —?, until a depth is reached where the velocity is double that at 
unit depth. The postulated depth-velocity function is then, of the type 


z=s(v—1)forilisv<y 


z= 1+ so(v — 01) — (k?/2)(v — 2)? for ») S v < 20. 


Let now —k? be so determined that crossing will be initiated at the depth 
corresponding to v = 22;. Substituting in Eq. (15) and performing the integra- 
tion, we obtain the condition 





So Tv Z 
= k*| — — sin $ |, 
[4(1 + 1/59)? — 1]/2 2 


0.95559 
[41 + 1/59)? — 1)? 
Thus for (a) so = 10 the critical value for k? is 4.86, and for (b) so =2, k? = 0.676. 


or 


2 

















12L 
Fig. 6, Critical ray paths corresponding to the velocity functions of Fig. 5. 


Larger values of k? in each case lead to a reversed segment in the travel-time 
curve. The velocity functions corresponding to these critical values are 
plotted in Fig. 5, and in Fig. 6 are shown the critical ray paths, which corre- 
spond. It is noteworthy that so gentle a variation in the velocity curves as 
is shown in Fig. 5 should be competent to produce crossing of the ray paths. 
If such gradually varying functions are able to produce this result we must antici- 
pate the frequent occurrence of the triple valued travel-time curve under conditions 
of seismic prospecting. 


43 The earlier computations of the heights reached by acoustic rays in the stratosphere dis- 
regarded, in effect, the triple valued nature of the travel-time curve. Thus Gutenberg (Gerland’s 
Beitrage Z. Geophysik 27, 22 (1930); 35, 46 (1932) assumed a zero range in the stratosphere for 
grazing angles of incidence, and, in essence, a “normal” travel-time curve of type I. This over- 
sight was pointed out by Whipple (Gerland’s Beitrage Z. Geophysik 31, 158 (1931)) who ob- 
tained somewhat higher altitudes for a given temperature on recomputing the data. Whipple’s 
method avoids consideration of the “reversed” travel-time curve by the use of the solution of 
the Abelian equation in the form of Eq. (8) instead of Eq. (9). 
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3.3 Interpretation of the triple-valued time-distance curve 


If the range, 4/2, and the time, ¢/2, to the vertex of a ray are written in 
the form 


Zn 
A/2 = f tan i dz 
0 


Zm 
/2 = dz/Vm sin i cos i 


then clearly 


Zm 
t = A/tm+ 2f (cot i/vm)dz 
0 


Um 
A/tm + 2 f (cot i/0m)(dz/dv)dv 
vo 

where for a given velocity function we may regard the last integral as a 
function only of vm. Hence 0A/dt=vm, or 5=vm. Thus this relation, which 
was previously obtained (Eq. (3)) from simple geometrical considerations 
holds in general, and in particular, for the present case of the triple-valued 
travel-time curve. (Also since v,, is not zero, a focus point in the distance will 
imply a focus point in the time.) The other formulae used in section I are 
also valid so that the previous analysis is essentially unchanged by the 
occurrence of a reversed segment. However, as has been said, for increasing 
Ym the reversed segment is described in the direction of decreasing A, so 
that its contribution to the integral in Eq. (9) is negative. 

A reversed segment, since it is associated with a velocity function increas- 
ing with depth, and is described in the negative sense, will necessarily be con- 
cave upwards. It meets its adjoining segments of the travel-time curve tan- 
gentially. 


3.4 A special case involving a discontinuous velocity function 


An instructive example of the interpretation of a triple-valued travel- 
time curve is afforded by the case of a contact between two uniform beds of 
different velocities. This case is very simple, and special methods for inter- 
preting its travel-time curve are in common use. Indeed, the fact that the 
general method of Herglotz-Wiechert is applicable to such discontinuous 
cases or even to the case involving the reverse segment, seems to have been 
generally overlooked.“ Therefore, despite the simpler special methods avail- 
able, it seems worth while to consider this case from the point of view of the 
general theory, especially since this discontinuous case may be considered as 


4 For example, Gutenberg (Lehrbuch der Geophysik, Borntraeger, 1929, p. 250) writes 
“The (Wiechert) method is valid only for cases in which the velocity does not change dis- 
continuously,” and similarly in his Hand. d. Geophysik, IV-1, p. 42, “the method breaks down 
when discontinuities of the first order exist, because then the travel-time curve consists of 
several disconnected portions.” Similarly, V. Conrad (Enzyklopidie d. Math. Wiss. 6, part 1, 
p. 449) writes “The Herglotz theorem is valid only when the rays . . . possess no focal surface.” 
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a limiting one representative of a family of similar continuous cases which are 
not as a rule amenable to treatment by special devices. 

The upper layer will be assumed to possess a thickness h, and a velocity 
Vo; the underlying material, a velocity Vi, Vi> Vo. The resulting travel-time 
curve, shown in Fig. 7, will consist of three segments, the initial straight line 
(1) of slope corresponding to the velocity Vo, the “reversed segment” (2), 
and the straight line (3) corresponding to the velocity V;. Theoretically, of 
course, the energy received along the straight segments is zero but to avoid 
such difficulties we may regard the example as the limiting case of one in 
which the velocity in the two regions is not constant, but increases infini- 
tesimally with depth so that a finite energy is received, which vanishes as the 








r 
 ] 
2 
7 
aL 3 
+ an 
+ V,= 1.5 Vo 
6 aN 
4- 
A 
2L- 1 
1 1 1 1 1 i l 4 J 
0 2 4 6 8 10 


Y, 


Fig. 7. Travel-time curve for the discontinuous “step” velocity function. 


travel-time curve approaches the straight lines shown. The reversed segment 
(2) in the limit coincides with the time-distance curve yielded by the wave re- 
flected at the surface z=h, and its equation is therefore: 


t = (A? + 4h?)!/2/y, 
and the apparent velocity, dA/dt, is 
BD = v9(4h?/A? + 1)!/2, (15a) 


The refracted travel-time curve begins at the point marked A, (Fig. 7) which 
corresponds to the angle of total reflection and therefore has coordinates: 


{ A = 2h/[(v1/00)? — 1]*/? 
t (2 hv/ v9) (,? —- V9") 1/2, 


II 


To the left of this point the curve corresponds to the reflected wave alone. 

The interpretation of this travel-time curve obtained in accordance with 
the general method is as follows. For any point on the segment (3), velocity 
V;, the corresponding depth is given by the formula (see Eq. 9) 
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| =) 1 a 1 
| TZ = f cosh~! — dA — f cosh~! ———_—_—__ wrod dA 
0 vo 2h/{(vy/r_)2—1Y2 vo(4h?/A? + 1)! 
where the second integral represents the contribution of the reverse segment. 
On substituting v;/m=k and k/(4h?/A*+1)!/?=x this becomes 
ah) (kP—1 ye * (cosh~' k — cosh x) 
= J cosh~! kdA +f ener see game 2hk*dx. (16) 
0 , (k? — 2) 3/2 
After an integration by parts this yields, as it should, 
{ wz = th. (17) 
) 
| The depth corresponding to any velocity on the reversed segment 1s also 


readily seen to be h. For at any point P on this segment, of epicentral distance 
A, and apparent velocity v,, the depth z, is given as before by 





Ap Y x v : 
T2p = f cosh7! va dA a J cosh =} — — cosh- 1 Vp t dA 
0 Vo Ap tT Vo v9(4h?/ A? + 1)1/2 | 
Put 
k’/(4k2/A? + 1)", 


ll 


Vp/%o = k’, and x 

Now by (Eq. 15a) 

v = v9(4h?/A? + 1)!/2 
| A = 2hvo/(v? — v9?)'/?. 
Hence 

Ap = 2h/(k”? — 1)*/2, 

Also A=A, implies 

x = k'(4h?/A,? + 1)'/= 1 
and A= implies 





x= k’. 


Thus 


2h (k?*—1)1/2 k’ (cosh 1 kb! oe cosh 1 x) Dhk’? 
TZ p = f cosh tl k’dA oe i ceed oS ED dx. 
0 ; (k” — x2)3/2 


This is precisely Eq. (16) above except that k’ replaces k, and the integration 
obviously yields the former result +z, = 7h, for all points P. Thus the depth 
interpretation for all points along the segments (3) and (2) is uniformly h. 

In summary, the interpretation yielded by this travel-time curve is the 
velocity function shown in Fig. 8. The depth h is yielded by the interpretation 
above, and the dotted portions follow from the fact that the straightness of 
segments (1) and (3) indicate the constancy of the velocity in the two regions 
involved. 


eee ee! 


3.5 Error involved in disregarding the travel-time curve at large epicentral 
distances 





For the discontinuous depth-velocity function of this example, the seg- 
ments (1) and (2) are asymptotic at infinity, and for functions of similar 
type, though not strictly discontinuous, the junction point of the segments 
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will usually occur at epicentral distances too large to be exhibited by the 
field data. The error involved in disregarding the contribution of these seg- 
ments at large epicentral distances is seen as follows. If Eq. (17) is written: 


TZ = thy, + [ 


where /, is the depth computed when the portions of the travel-time curve at 
epicentral distances greater than nh are omitted and C is the correction 
thereby involved, then from Eq. (16) 


’ 2hk?dx 
J (cosh—! k — cosh~! «)————__— 
k 


Al n24.1)"2} ( k2 mes x2) 3/2 


C 





k 
— (cosh k — cosh~! ——-—————— ) mn 
[(2/m)? + 1]? 


k* dé 
+h f ae : 
K/[4/n?+1) (kh? — §)*/3(— — 1)*/2 








O U 
; 
! 
| 
' 
\ 
' 
' 
ty 
Fig. 8. 


where the substitution =x? has been introduced in the last integral. Hence 


C = h |e — 2 tan (n/2k)(k? — 1 — 4/n?)1/2 


; k 
— n\ cosh! k — cosh=! ———_———_— } ]. (18) 
(4/ n2 + 1) 1/2 


From this formula the value of the relative error in the depth determination 
for the special case k=1.5 plotted in Fig. 7, and for several values of 7, is as 
follows: 


n 5 10 50 


(h— hy)/h 0.165 0.082 0.014 
Thus the error committed in neglecting the portion of the travel-time curve 
at distances in excess of ten times the thickness of the upper layer is small. 


3.6 Possibility of confusing curves of types II and III 


The rational interpretation of travel-time curves of type III may some- 
times be quite impossible in practice, due to the difficulty in identifying the 
reversed segment (2) and in constructing a correct travel-time curve from 
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the field data. The construction of the reversed segment involves the recog- 
nition of third arrivals of the wave pulses, but when the second arrivals corre- 
sponding to the extension of the (nearly) straight portions (1) and (3) of the 
travel-time curves can be recognized, the existence of the connecting seg- 
ment (2) is thereby established. It should then often be possible to establish 
the general form of the connecting segment with sufficient accuracy by the 
conditions that it must be concave upward and must meet the two known 
segments tangentially. 

In any event, it is of theoretical importance to distinguish between a 
travel-time curve having a discontinuity in slope at a focus point, as in Fig. 
3b, and one which has a reversed segment, as in Fig. 4, for the latter type of 
curve implies a much more rapidly increasing velocity with depth than the 
former. To illustrate this point suppose the later arrivals in the travel-time 
curve of Fig. 7 to be disregarded, so that the curve consists only of the 
(nearly) straight segments (1) and (3) up to their point of intersection. By 
the standard formula, the interpretation of this curve is 


z = (A,/z) cosh v/2, 
or 


v = v, cosh 12/A; 


where 4; is the abscissa of the point of intersection, and is therefore equal to 
2h|(v2/v; — 1) /(v2/21 +1) |*/2, where h is the depth yielded by the solution of the 
complete travel-time curve of Fig. 7. As mentioned in §2.2, at the point of 
intersection the apparent velocity must be assumed to change continuously 
from its upper value v2 to its lower one, 2. 

This interpretation, a hyperbolic cosine velocity-function, is quite dif- 
ferent from the step function yielded by the complete curve. Indeed, when 
v2=1.5 v;, as in the example of Fig. 7, the hyperbolic cosine formula yields 
z=3.04h, instead of z=h as previously obtained for the case of the complete 
curve, as the depth for which the velocity reaches the value v2. This factor of 
3 to 1 is especially noteworthy because the two travel-time curves considered 
are identical as to first arrivals, and the subsequent arrivals may easily escape 
notice on the record. The possibility of confusing these two types of curves 
and of obtaining mistaken interpretations is apparent. 


Part IV. DISCONTINUOUS TRAVEL-TIME CURVES CAUSED BY SLOWER 
SPEED BEDs 


A fourth type of time-distance curve is the discontinuous one which oc- 
curs when a slower speed bed is present. In this case, as illustrated by Fig. 9, 
the ray which enters the slower bed at grazing incidence is deflected down- 
wards and does not reverse until it has penetrated to such a depth that the 
velocity has again attained its previous maximum. This ray may emerge at 
an epicentral distance considerably beyond that of its neighbor, so that a 
dead zone in which no energy is received results at the surface. In such cases, 
as is well known, it is impossible to deduce from the travel-time curve unique 
values for the seismic velocities at depth. Nevertheless, the theory of this 
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type of curve supplies us with extremal solutions which may permit deduc- 
tions of satisfactory geological accuracy. In considering the theory of this 
case, the travel-time curve is as usual assumed to be completely known from 
the field data. The required degree of completeness in the definition of travel- 
time curves containing a gap is, however, by no means assurable from the 
experimental point of view. This difficulty, however, is properly disregarded 
in considering the narrower question of the theoretical possibilities in inter- 
preting the ideal travel-time curves themselves. 


4.1 Rearrangement of the velocity function 


In dealing with the slower speed bed, it is convenient to utilize the fol- 
lowing fact. The portion of a travel-time curve which is generated by rays 
which have penetrated below any reference level, 2, is unaffected by any 





| 
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Fig. 9a. and Fig. 9b. Discontinuous travel-time curve corresponding to slow speed bed. 





rearrangement or “shuffling” of the velocity function in the region above this 
level. This is perhaps obvious, but may be made clear as follows. The range 
x of any ray at its intersection with the level of reference 2, is expressible by 


the integral: 
8 
r= I tan i dz. 
0 


Here, tan i at any level is fixed along any ray only by the velocity at that 
level, through Snell’s law. Thus any rearrangement of the velocity function 
determines a corresponding rearrangement of the values of tan 7 and the 
area denoted by the integral remains unchanged by the process. Also, the 
angle of emergence just below the level z; is fixed for any ray by the velocity 
at this level. Hence for all such rays the initial conditions, i.e., the angle of 
incidence and the x coordinate, just below the level of reference are unaltered 
by any vertical rearrangement of the velocities above this level. Conse- 
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quently the geometry of the subsequent portions of the ray paths remains 
unchanged, as does the corresponding part of the travel-time curve. 


vw te 
4.2 Minimum thickness of slower speed horizon 


In general, the velocity will have arbitrary variations in the slow speed 
bed, subject of course to the restrictions that v< V where V is the previously 
attained maximum. But in accordance with the previous paragraph we may 
imagine the velocities in the slower speed bed rearranged to give a non- 
decreasing function, without altering the observed travel-time curve in any 
way. At the boundaries of the gap in the travel-time curve, the apparent 
velocities will be V, the seismic velocity at the boundaries of the slower bed. 
The velocity function in this bed is subject to the restriction that it yield the 
observed gap in the travel-time curve, which, let us say, has a discontinuity 
T in the time, and X in the distance. The two conditions imposed are, then 


T/2 = f dz/v cos i (19) 
21 


22 
f tan 7 dz (20) 
z1 


where 2; is the known depth to the topof the bed as deduced by interpreting 
the initial part of the travel-time curve to the left of the gap, and 2 is 
unknown. From these equations we may fix an upper possible limit for Ze, 
as follows: Eq. (19) may be written: 


T #2 (tani + coti X 72 coti 
a f ( at f ds. (21) 
2 Z V wd. ¥ 


x/2 











Instead of considering T fixed and z, unknown, let us consider for the moment 
Z2 to be constant, and seek the velocity function which minimizes T for fixed 
X. The above equation may be written: 


T X Vi- V)?\"/2 dz 
+f [1 — @/V)*]"? dz 








2 @ - v dy 


Euler’s test’ then provides the condition for a minimum 
[1 — (v/V)?]*/2/x = constant = (1/V) cot i. 


Hence for any value of X the minimizing ray path is a straight line, 7=con- 
stant. In view of the auxiliary condition (Eq. (20)) the necessary value for 


1 is: 
i = tan-! X/2(ze — 2;). 
If h is the thickness of the bed, the minimum time, 7°», is then, from Eq. (21) 
T mn/2 = (h/V)\tan i + cot i} 
(X/2V) [1 + 4h?/X?]. 
4 A. R. Forsyth, Calculus of Variations, Cambridge Univ. Press, 1927, p. 16. 


(22) 
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Thus for fixed X, 7,, and h increase or decrease together. Hence for a given 
discontinuity 7, the maximum thickness h is, as given by Eq. (22): 


j [—(+ x yy 1(X)1/2 (23 
1 = quent aeuies: om = 4(AL /3 . 
2\2 2V ; ) 


where L denotes V7'—X, and is the distance marked in Fig. 9. 
It follows that the velocity v,., which insures the maximum thickness for 
the slower layer has the constant value 








hes X/2 VX\"/2 
v% = Vsini=V ——— (=) (24) 
(XL/4 + X2/4)12 T 


4.3 Minimum thickness of slower bed 


The minimum thickness of a slower horizon consistent with observed dis- 
continuities at the gap in the travel-time curve is also easily found. In any 
actual problem, one may always fix a limiting value, v,, below which the 
velocity surely will not fall in any beds which can be present. This minimum 
velocity will evidently give a minimum thickness for the slower bed. How- 








Fig. 10. Ray path for minimum thickness of slower bed. 


ever, in order to satisfy the observed discontinuity, X, in the epicentral dis- 
tance, an infinitesimal thickness must be considered added to the lower part 
of the slow layer, in which the velocity increases very rapidly from the value 
v, to its value V at the termination of the slower bed. The ray path will then 
be as illustrated in Fig. 10, from which the discontinuity T is found to be: 





T 2hi X/2 — hy tani 
2 V sin 2: V 


Hence the minimum thickness, /;, is 
tan 7 L v1 


hy = ——{TV — X} = - 
2 2 (V2 — 9/2)12 





We have determined, then, the maximum and minimum thickness for the 
slower speed layer, and the values of the velocities which correspond. By 
using these results, interpretations for the remaining travel-time curve be- 
yond the gap may be obtained. 
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4.4 Interpretation of the travel-time curve beyond the discontinuity 


After interpreting the initial part of the travel-time curve up to the gap 
in the standard way, we may use this result to form two auxiliary depth- 
velocity functions. (1) The first function is a combination of this result with 
the solution corresponding to the maximum thickness of the slower speed bed, 
but rearranged or “shuffled” to give a non-decreasing function. (2) The sec- 
ond is the similar non-decreasing rearrangement formed for the minimum 
thickness and velocity of the slower bed. 

As has been said, the travel-time curve beyond the gap commences with 
a slope to correspond with the boundary velocity of the slower bed, V. Be- 
yond this, however, we may not in general say what its subsequent form or 
type will be, but for simplicity we will assume that no additional slower beds 
occur at depth, so that the subsequent travel-time curve will at all events 
correspond to velocities increasing in depth. 

We have, then, to consider the two known, but artificial, velocity func- 
tions (1) and (2) which are constructed as explained for the slower bed and 
the region above it; and the known second segment’ of the travel-time 
curve which refers to the region below the gap. 

From the known upper velocity functions (1) and (2) the corresponding 
travel-time curves may be computed. These curves will necessarily contain 
a “reversed” segment, asymptotic to a straight-line segment extending to 
infinity (as illustrated in section III) by virtue of the layer of constant ve- 
locity which enters the velocity-functions. In each case, of course, these 
travel-time curves will join the second segment of the original curve continu- 
ously for this condition was imposed in satisfying the discontinuities at the 
gap. On combining these curves, we obtain two continuous travel-time curves 
of the form shown in Fig. 11, of which the segments marked (2) are alike. 
These curves are now interpretable by the usual graphical method, and yield 
two extremal interpretations for the region below the slower bed, one corre- 
sponding to the assumption of maximum thickness, the other to the assump- 
tion of minimum thickness for it. 

Thus our complete interpretation of the discontinuous travel-time curve 
involves three parts; (1) the interpretation gained from the initial segment 
in the normal way; (2) a maximum and minimum interpretation gained from 
the gap, and (3) the corresponding two interpretations gained from the sec- 
ond segment as described above. The result is two velocity functions such as 
are sketched in Fig. 12. In practical cases it will probably often be found that 


16 This second segment, let is be noted, may be of a complex form, as produced, let us say, 
by a combination of segments of the normal and reversed type. (However, the occurrence 
no additional gaps or discontinuities has been postulated.) The possibilities regarding the form 
of this second part are in no way limited by the character of the known functions (1) and (2) 
as deduced for the region above the bottom of the slower bed. For in cases of velocities increas- 
ing with depth, as noted in paragraph 2.2, the derivative dy(A)/dd may be considered to vary at 
will from negative to positive values, and vice-versa, corresponding to crossing or non-crossing of 
the rays. Thus in such cases the type of travel-time curve traced by rays penetrating below 
a given level is essentially unrestricted by the velocity function above this level. 
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these two solutions are sufficiently alike so that their mean may be used with- 
out danger of an error significant from the geological point of view. 
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Fig. 11. Auxiliary travel-time curves for interpretation of region below slower bed. 


In case an additional slow speed bed is encountered at depth, the above 
procedure of course leads to four, instead of two, tentative solutions, etc. 
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Fig. 12. Extremal interpretations for case of slower bed. 


Part V. COMMENTS REGARDING THE APPLICATION OF THE ABOVE RESULTS 
TO A SPHERICAL EARTH 


With slight modifications the above discussion of seismic travel-time 
curves for a flat earth will in general also apply to the corresponding topics 
in the case of a spherically symmetrical earth. But an exception to the rule, 
and one modification of the foregoing results will be specifically mentioned. 
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5.1 Failure of the previous method in the case of a slower speed shell 


The analogous interpretations of the region inside a slower speed shell in 
a spherical earth may not be obtained by the method used in the case of a 
flat earth, i.e., by the device of rearranging the velocity function in a non- 
decreasing way. For in the spherical case the travel-time curve corresponding 
to the underlying region is altered by applying such a process to the region 
above. More specifically, the angular range traversed by a ray in passing be- 
tween the levels 7; and re is 


6 = -{ tan idr/r 
r1 


where tan 7 is given by the well-known relation, r sin i/v=1r9 sin i9/vp =con- 
stant along any ray. Here 1p is the radius of the earth, 7 the radial distance 
from its center, and the other symbols are as previously defined. Thus the 
expression for @ is of the form 


i oe, 





¢=— f(r, r/v)dr 
r1 
= f(r, o(r))dr. 
ry 
In general this is obviously not constant when ¢(r) is varied by the device of 


“shuffling” v(r); hence the former method fails in this case. 


5.2 Occurrence of the triple valued travel-time curve 





The criterion for the crossing of ray paths in a spherical earth which cor- 
responds to Eq. (16) above may be obtained as before by writing the half 
range to the bottom of a ray as an integral in the variable v, and differenti- 
ating with respect to the upper limit. Thus 


Am if tan idr 
2ro . r 


ft (%m/Um)(dr/r) _ 
ro (4/0)? — (4m/m)?]!? 
f Ym v(dr/dv)dv 


Yo [(omt/ 1m)? ne y?}!/2 














and 











1 oa =| (dr/dv),=r, iZ (d*r/dv*)do | 
| (ro0m/Tm)? = v2] v6 [(ro0m/Tm)? om v?| 1/2 , 


Thus the condition for the occurrence of the triple valued time-distance 
curve is: 


2ro dim Tm? 








Ir/ dv) pa» Ym 1*7r/ dv?) dv 
(dr/dv) rar, f (d?r/dv*) wis (25) 


[(ro0m/1m)? _ Yo? ]1/2 m [(ro0m/Tm)? = y? ]1/2 


and vice-versa. 
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The application of this criterion to a velocity function which approxi- 
mates that for the earth seems to indicate that there is‘a much lesser likelihood 
for the occurrence of a triple valued curve than was noted above in the case 
of seismic prospecting. For example, let the seismic velocity function in the 
earth be taken to be that given by Gutenberg.'? Between the surface and a 
depth of 1200 km a linear variation, v = 32.7 —24 r/ro prevails, corresponding 
to a velocity of 7.9 km/sec. at 60 km depth and 12.2 km/sec. at 1200 km 
depth. 

Now consider a slightly different function commencing with the same in- 
itial velocity, having a constant slope 10 percent less than the above until 
the depth corresponding to v=10 is attained, and thereafter possessing a 
constant second derivative d’r/dv? until the velocity v=12 is attained. Fur- 
thermore, let the value of d?r/dv? =k? be so chosen that the “reversed” travel- 
time curve is incipient when v=12. The postulated velocity function is, then 


Vv 


29.3 — 21.6r/ro when v < 10 
and, 


r/ Vo 


when v>v. Here v is taken as 10. From Eq. (25), the condition on k?/rp is, 


0.892 — (v — )/21.6 + (k?/2r9)(v — 21)? (26) 








1 k? Tm 10 1m 
0. =—- ——— + —) sin! — — ging? — — 
21.6[(12r0/rm)? — (7.67)?]!/2 or ro 12 ro 


where rp and k?/r also satisfy Eq. (26) formed for v=v,,. By trial, the solution 
of these equations is found to be k?/r9=0.0175. The resulting velocity func- 
tion is 


r/ro 1.0 0.893 0.855 0.835 0.833 


v 7.67 10 11 12 13 





This is plotted in Fig. 13 in comparison with Gutenberg’s function. There is 
also shown a third function which is obtained as just described, with the dif- 
ference that the junction of the straight and curved portions is supposed to 
occur at v=9 instead of v=10. In both cases, the condition is imposed that 
incipient crossing occurs for rays which bottom at the depth where v=12. 
These altered velocity functions differ sufficiently from the actual one as 
given by Gutenberg so that it is unlikely that this requisite difference can 
really exist. The conclusion thus seems justified that in general the travel- 
time curves will not be triple valued from such a cause, or if so, at least not 
to such an extent as to affect significantly their interpretation at depth. This 


of course is in agreement with the experimental evidence and the usual con- 
clusions relative to this point. 


CONCLUSIONS 


It appears that the interpretation of seismic travel-time curves in horizon- 
tal structures has a satisfactory theoretical basis, under the assumptions 


17 B, Gutenberg, Der Aufbau der Erde, 1925, p. 31. 
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made. The classical method of Herglotz-Wiechert may be extended to the 
general case, involving multiple valued travel-time curves, discontinuities in 
the velocity depth function, and also, in a certain sense, to cases involving 
slower speed beds. The experimental definition of the travel-time data is not, 
however, always sufficiently complete to fulfill the needs of the theory. The 
case of the slower speed bed, and some cases in which an accelerated increase 
of the velocity occurs with increasing depth, offer experimental difficulties in 






Z_ Actual Function 


Assumption a per Gutenberg 


~—Assumption m4 





6 | er 1 } 
|-O 09 08 O7 
Tre 








Fig. 13. Example showing velocity functions which produce triple 
valued travel-time curves in case of a spherical earth. 


the complete identification of the travel-time curve. The errors introduced 
by deficiencies in the experimental data will in certain instances cause the 
velocity to appear to vary more gradually than is really the case. 

The writer wishes to express his appreciation of the kindness of Professor 
P. Franklin in reading the manuscript. 
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The Measurement of the Phase Angles of Shielded Resistors 


By Leo J. BERBERICH* 
Electrical Engineering Department, Johns Hopkins, Baltimore, Md. 


(Received October 3, 1932) 


In this paper is discussed the problem of measurement of the small, but not al- 
ways negligible, phase angles of resistors, with special reference to shielded resistors 
of 1000 ohms and above. Several types of shielded resistors are described briefly, 
among which is the ideally shielded, low-phase-angle resistor, recently constructed 
by the writer. The substitution method developed is one in which the resistor, the 
phase angle of which is desired, is compared, by means of bridge circuit, toa reference 
standard of a simple geometric configuration and of approximately the same re- 
sistance as that of the resistor. The results of measurements, at frequencies of 500 
and 60 cycles, on three ideally shielded resistors are given. By a series connection of 
resistors it is possible to extend the range of the method to 50,000 ohms. 


1. INTRODUCTION 


S OUR standards of accuracy progress precision of measurement be- 
comes increasingly more important. In precise alternating current meas- 
urements cases often arise where the small phase angles of resistances cannot 
be neglected. Because of the generally conceded difficulty of the measure- 
ment of such phase angles only meager information appears in the literature 
concerning this problem. Noteworthy is the work of Grover and Curtis! of 
the Bureau of Standards, and as a general treatment of the subject it may be 
considered monumental. Of the early work done in this connection that of 
E. Giebe’ and that of Wien and Prerauer* deserve consideration. In both the 
former and the latter a modification of Maxwell’s method for the comparison 
of two small inductances was used. Most of the work done in this field up to 
the present time was carried out at frequencies in the audible range and only 
very little at frequencies as low as 60 cycles per second. 

The question of choosing a method for this type of measurement is not 
settled by taking into account merely the magnitudes of the residual induc- 
tance and capacitance. The value of the resistance with which they are as- 
sociated must also be considered. Thus the “time constant,” which is the ratio 
of the “effective inductance” (resultant effect of residual capacitance and in- 
ductance) to the resistance, is a determining factor as well as the magnitude 
of the resistance. The problem of measuring the phase angle of a low re- 
sistance having a given time constant is, in general, different from that of a 
high resistance having the same time constant. In low resistance coils, 1 to 


* Researsh and Development Department Vacuun Oil Company, Inc., Paulsboro, N. J. 

1F, W. Grover and H. L. Curtis, The Measurement of the Inductance of Resistance Coils, 
Bull. Bur. Stds. 8, 463 (1913). 

*E. Giebe, Prazisionsmessungen an Selbsinduktionsnormalen, Zeits. f. Inst. 31, 6-20, 
33-52 (1911). Also see B. Hague, A. C. Bridge Methods, 2nd Edition, p. 216. 

3 O. Prerauer, Ann. d. Physik 53, 772-774 (1894). 
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10 ohms, the magnetic field close to the wires is not entirely eliminated and 
gives rise to a small residual inductance. There is also, however, an electro- 
static capacitance between adjacent turns and adjacent layers which tends 
to neutralize the inductive effect of the magnetic field. In 100 ohm coils, bi- 
filarly wound, it has been found for low frequencies that these two effects 
approximately balance each other, and thus they are essentially nonreactive. 
When we come to resistance of 1000 ohms or above, however, we find that the 
effect of electrostatic capacitance predominates. If we consider this distrib- 
uted shunt capacitance as lumped and connected across the terminals of the 
resistance, it is a simple matter to show that the capacitance effect is propor- 
tional to the product of the square of the resistance R, and the capacitance 
C which shunts R. Hence, the effect of capacitance increases as the square of 
the resistance and the higher resistances can be expected to be capacitative 
rather than inductive. 

Grover and Curtis divided their investigation of this problem into two 
parts; the first part treats methods for measurement of low resistance coils, 
and the second part, methods adapted to coils of resistances greater than a 
1000 ohms. The reason for this division is quite logical. In measurements on 
low resistance coils the inductance of the leads and all mutual inductance 
effects must be considered or else shown to be negligible. For coils of high 
resistance the inductance effects become negligible, and the capacitance 
effects between the various parts of the measuring circuit and to earth be- 
come important. The two cases thus differ essentially in that the first involves 
the measurement of an extremely small inductance and the second, that of 
an extremely small capacitance. In this paper, however, only the second 
case will be treated since we are here concerned mainly with shielded resist- 
ances of 1000 ohms and upward. One of the methods used by Grover and 
Curtis is modified by introducing certain novel features and simplifications. 

In order to simplify subsequent computations the residual capacitance 
will be regarded as a negative inductance insofar as its effect on the phase 
angle is concerned. It is easy to show that the tangent of the phase angle of 
a coil of resistance R, inductance L, and shunt capacitance C is approxi- 
mately the ratio of 2xf (L — CR?) to R. Now it is convenient to call (L — CR?) 
the effective inductance, L’, of the coil. L’ will be either positive or negative 
according as the effect of inductance or capacitance predominates. Hence, 
the effective inductance may be defined as, the inductance which when multi- 
plied by 27f and divided by R gives the tangent of the observed phase angle. 
The phase angle is positive when L’ is positive and negative when L’ is nega- 
tive. 

2. SEVERAL TYPEs OF SHIELDED RESISTORS 


In general every resistance contains three parasitic properties; namely, 
self-inductance, self-capacitance (capacitance between the various parts of 
the coil), and earth capacitance. Of these three, the last named is the most 
troublesome. The first two can be made small by the proper choice of a wind- 
ing for the particular coil under consideration.‘ It is, however, not quite such 
a simple matter to reduce the effect of earth capacitances. Such capacitances 
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will cause the phase angle to vary with the position of the coil with respect 
to earth and also with the potential of nearby objects in the laboratory. The 
only remedial measure which is at our disposal is to completely shield the 
coil in a metal case maintained at a definite potential. This may be done in 
a number of ways. 

The simplest method of shielding is illustrated in Fig. 1A. Here the re- 
sistance is completely surrounded by a metal container which is connected 
to one end of the resistance in order to fix its potential. This, however, has 
the effect of increasing the phase angle because of the close proximity of the 
shield to the resistance. The distributed capacitances from the resistance to 
the shield are in general greater than similar capacitances to earth would be 
without the shielding. But such a shield will fix the phase angle and once it 
is measured it can be allowed for in the computations. 

As the resistance within the shield increases the effective inductance rises 
rapidly, approximately as the square of the resistance. The effective induc- 
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Fig. 1. Three types of shielding. 


tance can be decreased markedly by subdividing the resistance into sections 
and surrounding each section with a shield. Fig. 1B shows the “working re- 
sistance,” R,,, divided into two sections. The shields are maintained at a 
definite potential by connecting them to an auxiliary circuit R,. This auxili- 
ary circuit is known as the “guard circuit.” Shielded resistors of this type have 
been described and mathematically analyzed by Orlich,5 Kouwenhoven,‘ 
Hiecke’ and Silsbee.* The mathematical analysis of a series connection of a 
number of these resistors shows that the following conditions must be fulfilled 
in order to obtain a minimum phase angle for the combination: (1) the po- 
tential of the shield of every unit except the first and the last must be equal 
to that of the midpoint of the resistance within it, (2) the potential of the 
shield of the first unit (low-voltage end) must be below the potential of the 
midpoint of the resistance within it by an amount equal to one-twelfth the 
voltage across the terminals of this unit, and (3) the potential of the shield 
of the last unit (high-voltage end) must be above the potential of the mid- 
point of the resistance within it by an amount also equal to one-twelfth of 


*H. L. Curtis and F. W. Grover, Resistance Coils for Alternating Current Work, Bull. Bur. 
Stds. 8, 495 (1913). 


5 E. Orlich and H. Schultze, Archiv. f. Elek. 1, 1 (1912). 
6 W. B. Kouwenhoven, Washington University Studies 1, part 1, No. 2, p. 143 (1914). 
7 R. Hiecke, Electrotech. u. Masch. 33, 505 (1915). 


8 F. B. Silsbee, A Shielded Resistor for Voltage Transformer Testing, Bur. Stds. Sci. Paper 
No. 516, (1925). 
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the voltage across the terminals of this unit. In order to be able to interchange 
the various units of the combination or to tap off a portion of the total re- 
sistance it is often convenient not to fulfill conditions (2) and (3) and to ad- 
just the potential of every shield equal to that of the midpoint of the 
resistance within it. This is done in the case of the shielded resistor in use at 
the Bureau of Standards.* The result is a slightly higher phase angle than 
is the case when all the conditions are fulfilled. Thus, if the circuit is so con- 
structed that R, is equal to R, for every unit and the shield of every unit 
connected to the midpoint of R,, we have the second case. It is obvious that 
when this is done the capacitance currents from R, to the shielding are bal- 
anced; that is, just as much current will flow away from R,, to the shield in 
one-half of the unit as will flow back to R, in the other half. There will be 
no loss of current, theoretically at least, from R, to R,. This type of shielding 
results in a marked reduction in phase angle, yet it is not negligible and in 
many cases must be measured. 

If we carry the subdivision of the resistor to the limit where each element 
of the resistance R,, is shielded by a corresponding element of R, we have 
what is known as the ideal type of shielding. Fig. 1C indicates schematically 
such a resistor. The guard circuit should be constructed so as to be coaxial 
with, and completely surround the working circuit. Also, the resistance must 
be uniformly distributed and the same for both circuits. When this is the case, 
even though the distributed capacitances between R,, and R, exist, there will 
be no voltage across them because the potential is the same along correspond- 
ing points of the two circuits. A special case, one end grounded, of this ideally 
shielded resistor has been mathematically analyzed by Silsbee* who showed 
the phase angle to be extremely small, though not zero, because of the dis- 
tributed earth capacitances from the guard circuit to ground. A number of 
such ideally shielded resistors have been actually constructed by the writer 
and will be described in a later paper.* The ultimate aim in the development 
of the method to be described was to obtain a measurement at 60 cycles of 
the extremely low values of effective inductance which one may expect with 
this type of shielding. The results of a series of measurements on a 10,000 ohm 
unit, 5000 ohm unit, and 1000 ohm unit will be given in this paper. 


3. GENERAL PRINCIPLES 

Practically any method used to measure the effective inductance of re- 
sistance coils resolves itself into a substitution method. This affords the sim- 
plest method of eliminating the effect of parasitic influences in the circuit. 
Other methods have been used with some degree of satisfaction, but in the 
writer’s opinion the substitution method is perhaps the easiest to use, does 
not require any very special apparatus, and gives the most accurate results. 
By a substitution method is meant one in which a reference standard of cal- 
culable residuals and of approximately the same resistance as the unknown 
is inserted in one arm of a bridge network and a balance obtained; then the 


® W. B. Kouwenhoven and Leo J. Berberich, “A Standard of Low Power Factor.” To be 
presented at Winter Convention of A.I.E.E., New York, N. Y., (1933). 
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standard is replaced by the unknown and the network again adjusted to a 
balance by varying one arm only, all other arms being kept fixed during the 
substitution process. The amount of the change required in this one arm is 
usually so small, provided the resistance of the standard is adjusted to very 
nearly that of the unknown, that the change in the parasitic influences in the 
arm is negligible. The parasitic admittances of the fixed arms are automat- 
ically eliminated since they enter into the computations only as ratios. 


Standards of reference 


The standards used in this work were constructed in the form of two paral- 
lel wires of high resistance material. For 5000 ohm and 10,000 ohm standards 
special Nichrome wire was used, having a diameter of approximately 1.4 mils 
(0.035 mm) and having a resistance of approximately 250 ohms per foot. 
The spacing of the wires was maintained by means of four small brass pins 
embedded in two hard rubber spacing blocks, the latter serving as supports 
at the ends of the standard. The two parallel wires were joined at the far 
end and the open end was placed near the bridge circuit. At the open end, the 
standard was fitted with short copper leads which dipped into mercury wells 
placed as closely as possible to the brass pins in the spacing block. The leads 
from the bridge circuit made contact with two other mercury wells, one be- 
side each of the two for the standard. This permitted the insertion of the 
standard in the bridge by merely dropping two small copper links, one each 
into the two pairs of wells. The two spacing blocks were placed at such a dis- 
tance that the wires were supported without appreciable sagging. The 10,000 
ohm standards were approximately 20 feet long. In order to obtain a standard 
of more than a few feet in length for the 1000 ohm case, No. 40 A.W.G. man- 
ganin wire was used instead of the high resistance Nichrome. It is necessary 
that the standard have a considerable length in order to make the end effects 
negligible. 


Computation of residuals for standards 


The residuals of a standard such as described above can be easily com- 
puted. The inductance due to the magnetic field of the two parallel wires, 
having a radius r, length of a single wire /, and spacing d, is given with suffi- 
cient accuracy by the formula!® 


L = 4l(log. d/r + }) X 10~° henries 


provided / is greater than 10d and d is greater than 10r. 

The capacitance between the wires is not of great importance in standards 
of low resistance, but in the case of high resistance such as used here, it be- 
comes the predominating factor. The capacitance can be computed only when 
certain conditions are fulfilled which will be examined now. In general, when 
it is desired to compute accurately the capacitance between two conductors 
we must take into account their capacitances to earth or devise a means of 
eliminating their effects. 


10 Russel, Alternating Currents, Vol. 1. 
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Fig. 2 shows two wires in section suspended at equal distances from the 
earth and indicates the three capacitances which must be considered. C2 is 
the capacitance between conductors 1 and 2, and C; and C2 are the earth ca- 
pacitances. Grover and Curtis' have shown that for an important special case, 
the capacitance of a system such as shown in Fig. 2 is 


C= Cie +C,/2. 


This special case requires, (1) that C; be equal to C2, or that the two wires 
be of equal diameters and at equal distances from the earth; and (2) that V; 
be equal to — V2, where V; and V2 are the potentials of conductors 1 and 2, 
respectively. The quantity C is the capacitance in the usual sense, between 
the two wires and is equal to the joint capacitance of Cz in parallel with C; 
and C: in series. Orlich, who first made the analysis, gives C the name “Be- 
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Fig. 2. A parallel wire standard in section showing the various 
capacitances involved. 


triebskapazitit” or “working capacitance.” This capacitance obtains only 
when the condition is fulfilled that the potential of one wire is just as much 
above ground as the other is below, or Vi= — V2. 
The working capacitance of the two parallel wires at equal distances from 
earth is!® 
C Kl X 10-* 
3.6(log. (d/r) — d?/4h? + d*/32h4*+.---) 


where K is the dielectric constant of the medium; / the length of each wire; 
r the radius of the wires; d the distance between centers of wires; and h the 
height of wires above ground. The correction terms in d/h are usually neg- 
ligible provided h>10 d. 

The contribution to the effective inductance of the standard by the work- 
ing capacitance C is not —CR? but —CR?/3 approximately. This is occa- 
sioned by the fact that one pair of ends of the parallel wire standard are 
joined together, and therefore the voltage between the wires varies from zero 
at the closed end to the total voltage across the standard at the open end. 
C is the capacitance which would be measured if the two wires were discon- 
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nected at the far end. The rigorous expression for the impedance of a stand- 
ard with one pair of ends joined is!! 


R + iwl\'/? 
Z= (<=) tanh [iwC(R + iwL) ]}/2. 


1 


Since both C and L are small this may be expanded into a rapidly converging 
power series 


Z = (R + iwL)|1 — (})iwC(R + iwL) + (2/15)iwC(R + iwl) +--- | 


where 7 is the complex operator and w is 27Xfrequency. Dropping terms 
involving the product of C and Z or their squares, this becomes 


Z=R+ iw(L — CR*/3). 
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Fig. 3. The circuit for measuring the phase angle of an ideally 
shielded resistor. 


Thus, it is seen that the effective inductance of a standard as used here is 
L’ = (L — CR?/3) 
where L and C are computed by use of the formulas given above. 


4. DESCRIPTION OF THE BRIDGE CIRCUIT 


A schematic diagram of the circuit used in the measurement of the phase 
angle of the ten thousand ohm ideally shielded resistor is shown in Fig. 3. 


" B. Hague, A.C. Bridge Methods, 2nd Edition, p. 37. 
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The circuit used for the units of lower resistance was essentially the same with 
the exception that the equal arm bridge was departed from. The basic ele- 
ments of the bridge reduce to an ordinary four arm network. The other arms 
are added merely to maintain the shielding at the proper potentials and to 
furnish a convenient ground for the circuit. As previously mentioned this is a 
substitution method where the substitution is carried out in arm A. 

Referring to Fig. 3, R, represents the working coil of the unit whose phase 
angle is to be measured. R,, can be inserted into or removed from the circuit by 
means of the links aa’ which fit into the mercury wells. The shielding resistance 
R,, which completely surrounds R,, is always in the circuit and is connected 
in series with a variable resistance R,’. R, and R,’ in series are connected di- 
rectly across the source. The parallel wire standard, Rs, is so arranged that 
it can also be inserted in or removed from arm A by the use of links bb’. Arm 
B contains a fixed shielded resistor having an ohmic value of 10,000; the 
shield is connected to one side of the source. Arm D consists of another 10,000 
ohm shielded resistor which has connected in parallel with it a variable air 
capacitor. Both the shield of the resistor and that for the capacitor are con- 
nected to the other side of the source. Arm C contains a variable shielded 
resistor of the decade type which is also shunted by a variable air capacitor 
with maximum capacitance of 250 mmf. This resistance box is equipped with 
six decades of resistance one of which permits a variation of resistance as 
small as 0.01 ohm, the total being 10,000 ohms. 

W, and We, which are commonly known as the Wagner arm, were made 
up of three tube rheostats adjusted by means of a Wheatstone bridge to one- 
tenth the resistance of the arms B and D across which this arm is shunted. 
The purpose of the Wagner arm is to maintain the detector shielding at a 
potential equal to that of the detector lead within it and also makes possible 
the grounding of the bridge at the proper point. In order to make the phase 
relations of W2 the same as those of arm D of the bridge it should be shunted 
by a capacitor. It has been found both experimentally and mathematically 
by W. B. Kouwenhoven and A. Bands!” that a small unbalance as to the 
phase of the shielding potentials produces a negligible effect on the main 
bridge, provided the insulation resistance between the shield and the lead it 
contains is infinite. Furthermore, even if the potential of the detector shield 
should differ a small amount in phase from that of the detector lead no error 
is introduced since arms B and D remain fixed during the substitution process 
as well as the Wagner arm. Thus, the omitting of a capacitor across W2 is 
doubly justified. 

W, is divided into two parts and at the junction of these the bridge is 
grounded. This in effect places the midpoint of the parallel wire standard, R,, 
at ground potential without actually grounding the standard itself. Thus, the 
condition is fulfilled that one wire of the standard be just as much above 
ground potential as the other is below it, which permits the calculation of 
the capacitance of the standard. This brings about a balance in the capaci- 


12 W. B. Kouwenhoven and A. Bands, Jr., A High Sensitivity Power Factor Bridge, A.1.E,.E 
Trans. 51, 202 (1932). 
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tance currents to earth just as was the case in the shielded resistor illustrated 
by Fig. 1B. 

In order to facilitate the adjustment of the currents in the shielding coil 
and the working coil of the resistors, provisions were made so that a milli- 
ammeter could be inserted in series with either by the use of a simple switch- 
ing arrangement. The current in R, was adjusted by varying the voltage on 
the bridge and that in R, by varying R.’. 

The alternating e.m.f. was supplied to the bridge through a two-to-one 
step-up transformer. The midpoint of the primary was grounded to eliminate 
capacitance effects from the source. Measurements were made at two fre- 
quencies, namely 500 and 60 cycles. A small aeroplane generator served as a 
500 cycle source, and a regular 60 cycle alternator, driven by a direct current 
motor supplied from a storage battery, furnished the 60 cycle e.m.f. At 500 
cycles a Baldwin telephone receiver unit of approximately 1700 ohms resist- 
ance was used as the detector. The telephone unit was completely shielded 
by means of a metal container which was connected to the shielding of the 
detector leads. The effects of earth capacitance to the observer were avoided 
by listening at the far end of a piece of rubber tube, the other end of which 
was connected to the diaphragm chamber. For the 60 cycle measurements a 
high resistance (850 ohms), high current sensitivity (0.25 microamperes per 
cm at a meter scale distance) vibration galvanometer served as a detector. 
This instrument was recently developed by the Leeds and Northrup Com- 
pany. The balance was detected by means of a telescope placed five meters 
from the galvanometer. 

Measurements were made on three of these ideally shielded resistors, 
ranging from 1000 to 10,000 ohms. Fig. 3 gives the circuit for the 10,000 ohm 
unit where the resistance in each arm of the main bridge is approximately 
10,000 ohms. For measurements on the other two units the circuit was modi- 
fied only slightly. Arms C and D were always left the same as well as W2 of 
the Wagner arm. The appropriate standard of very nearly the same resistance 
as that of the unknown, together with the unknown, were placed in arm A so 
that one could be substituted for the other. The resistance in arm B was made 
approximately equal to that of the unknown. The resistance of W; was ad- 
justed to one-tenth that of arm B. The ground was always placed at the mid- 
point of W; to bring the midpoint of the standard to ground potential. 

In making the measurements the bridge arms B and D (Fig. 3) remained 
fixed throughout the substitution process. The capacitor in arm D was used 
merely as a ballast to bring the reading of that in arm C to a convenient part 
of the scale. The unknown and the standard were alternately inserted in arm 
A, and a balance was obtained each time by adjusting the resistance and 
capacitance of arm C only. If the resistance of the standard is adjusted closely 
to that of the unknown the required change in resistance of arm C during the 
substitution process is so small that the change in the residual effective in- 
ductance of this resistance is negligible. No particular care need be exercised 
in adjusting the resistance of the standard to the same value as that of the 
unknown when the effective inductance of the resistance of arm C is known 
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for every setting by a previous measurement. This phase angle can easily be 
taken into account in the computations as will be shown. 


5. THEORY OF THE BRIDGE CIRCUIT 
In the solution of the bridge equation the following notation is used: R, 
is the resistance of parallel wire standard; L, the effective inductance of paral- 
lel wire standard; R, the resistance of working coil of unknown; L,, the ef- 
fective inductance of unknown. R, and C, are the resistance and capacitance 
respectively, required in arm C to obtain a balance with the standard in cir- 
cuit. R, and C, are the resistance and capacitance respectively, required in 
arm C to obtain a balance with the unknown in circuit. L, is the effective in- 
ductance corresponding to setting R, in arm C. L, is the effective inductance 
corresponding to setting R, in arm C;i the complex operator; Z,= Ry +iwLw, 
or R,+iwL,; Z,=Ro; 
‘ Ro R, 
Z.=— ——— GF see 
1+ wR, 1 + 1.C.R, 
r Ra 
aa 
1 + tw aRa 


Now for zero current through the detector, we have 





Z./la = Za/Z» = constant (1) 


where the subscripts refer to the various arms of the bridge of Fig. 3. Eq. (1) 
is equated to a constant because Z, and Z, remain unaltered during the meas- 
urement. Hence, we can equate the ratio of Z, to Z, with the standard in the 
circuit to the same ratio with the unknown in the circuit. Doing this we have 


Z. R, 1 R 1 
a. (2) 


Le 1+ wC,R, R, + wl, 1+ wR, Rw + wl 
Solving for L,,, we obtain the formula which is used to compute the effective 
inductance of the unknown in terms of the known constants of the circuit, 


Ly = RAL/R. + RC. — RuCz), (3) 


a result sensibly independent of the frequency. 

In the above analysis it was tacitly assumed that the resistance in arm C 
has no phase angle which is not true. But, since this phase angle enters as 
that of the difference of two settings which are nearly equal when the resist- 
ance of the standard is adjusted closely to that of the unknown, Eq. (3) is 
true for this special case. If, however, the phase angle of the resistance of 
arm C is known, Eq. (3) can be made general by adding two correction terms. 
In the latter case it is not essential to adjust R, closely equal to R,,. With the 
addition of the correction terms Eq. (3) becomes, 


Ly = Ro[L./Re + (RCo + Le/Re) — (RuCu + Li /R.) | (4) 


where L, and L, must be obtained from previous measurements on the re- 
ststance in C. 


a 
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6. PRECAUTIONS OBSERVED IN MEASUREMENTS 


In the measurement of such minute quantities as are involved here, ex- 
treme care must be exercised to avoid the introduction of errors comparable 
in magnitude with the quantity to be measured. A small movement of any of 
the leads in the bridge proper is sufficient to give rise to a large error. To 
guard against this all leads were rigidly fastened. The effect of capacitance 
to the observer from the detector or the main bridge arms is sufficient to give 
rise to a similar error. This required careful shielding of the detector and 
bridge arms as well as remoteness of the observer from the detector. With 
proper shielding a movement of the observer merely changes the capacitances 
of the shields to earth and its only effect is a change in the load on the source. 
When the telephone detector already described was used listening was done 
at the end of a five foot rubber tube. The 60 cycle detector was at a distance 
of more than five meters from the observer and balance was detected by the 
use of a telescope. 

The Wagner arm, W,; and We2, may be another troublesome cause of 
error if improperly placed. During the course of the development of this 
method, the Wagner arm and source were connected to the points P and Q 
and the detector to points M and JN instead of as shown in Fig. 3. Hence, the 
Wagner arm, instead of being across the fixed arms B and D, was across two 
arms, one of which is varied in the substitution process. When C is changed, 
the Wagner arm remaining fixed, it is evident that the potential of the de- 
tector lead is shifted with respect to the potential of its shield. Also, a few 
ohms change in C will shift the point on the bridge which is at ground poten- 
tial. These two disturbances were responsible for inconsistent results which 
at first were found difficult to account for. Therefore, it may be said that the 
bridge will operate satisfactorily only when the Wagner arm is in the position 
indicated in Fig. 3. 

The success of this method depends on the use of a relatively high voltage 
as a source. A high voltage was especially necessary at 60 cycles in order to 
obtain the required sensitivity. The limiting factor in raising the voltage was 
the current carrying capacity of the parallel wire standard. The voltage could 
not be carried above 300 volts without obtaining annoying fluctuations due 
to heating in the standard. 


Check on accuracy 


In order to check the overall accuracy of the set-up two parallel wire 
standards of the same resistance but with different spacings were constructed. 
The spacing between wires of one of these was 1 cm and that for the other 
was 20 cm. This may be said to represent two extremes of spacing. The ef- 
fective inductance of an arbitrarily chosen shielded resistor of considerably 
higher phase angle than may be expected of the ideally shielded resistor, was 
measured first by using one of these two standards and then the other. The re- 
sistor was of the type illustrated in Fig. 1B and possessed an effective in- 
ductance of approximately —0.9 mh for a 10,000 ohm unit. The circuit of 
Fig. 3 was adapted to this resistor by replacing R, with it and omitting R, 
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and R,’. The shield was connected to the ground point of the Wagner arm 
which fixes the potential of the shield and makes it equal to that of the mid- 
point of the resistance within it. In this preliminary test the effective in- 
ductance as obtained from the two sets of data differed by less than ten 
percent, or a time constant of less than 10-* seconds. This may be looked 
upon as a good agreement, considering the severity of the test and the small- 
ness of the time constant (9X 10~-* sec.). A large portion of the discrepancy 
in the results obtained with the two standards was probably due to errors in- 
volved in the measurements of the physical dimensions of the standards. 

For the purpose of gaining a more definite idea as to the magnitude of 
such errors, liberal uncertainities were assumed in the measurement of dimen- 
sions of the standard. Two complete computations of the effective inductance 
were carried through with such values as would give the extreme difference 
in the final result. For the one centimeter spacing the discrepancy between 
the two extremes was found to be about four percent of the total effective 
inductance of the standard. In a similar computation for the 20 cm spacing 
it was found to be only one percent. The difference in the two cases is mainly 
due to the accuracy with which the spacing can be measured. It is reasonable 
to assume that a 20 cm spacing can be measured with greater relative accu- 
racy than a 1 cm spacing. These computations, therefore, indicate the im- 
portance of measuring the physical dimensions of the standards accurately 
and especially the spacing between wires. 


7. MEASURING PROCEDURE 


Since the ideally shielded resistors were constructed so that the resistance 
of the shielding coil is equal to that of the working coil, a minimum of ef- 
fective inductance may be expected when the currents through the two coils 
are equal. Actually, this is not quite the case because of end-effects, and the 
always existing self-inductance and self-capacitance. In order to gain some 
definite information as to the nature of the variation of the effective induc- 
tance of the working coil with variation of the current through the shielding 
coil, a series of measurements was made in which the current through the 
shielding coil was changed from zero to approximately twice the fixed current 
through the working coil. The current through the working coil is determined 
by the maximum permissible voltage on the bridge commensurate with heat- 
ing. The current through the shielding coil was adjusted to any desired value 
by varying the resistance R,’ which possessed a range of from 1 to 100,000 
ohms. 

It might be mentioned here that the resistance of the milliammeter used 
was of the order of several hundred ohms, and thus small compared to the 
resistance of the bridge arms with which it was connected in series to measure 
the current. This method of measuring the current was accurate enough for 
the purpose in hand. 

Again referring to Fig. 3, with R,, in circuit and the standard disconnected, 
the current through R,, was brought up to a fixed value by increasing the 
voltage on the bridge. Then the current through R, is adjusted to the desired 
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value by varying R,.’. This being done, the milliammeter is taken out of the 
circuit, and a balance is then obtained by use of the resistance and capaci- 
tance in arm C. It may be found necessary to adjust the capacitor in arm D 
to bring the balance point of the capacitor in arm C to a convenient part of 
the scale. Once the latter adjustment is made arms B and D are left un- 
touched. This balance results in a resistance reading R, and a capacitance 
reading C,, in arm C according to the notation used. Now by the use of links 
aa’ and bb’, R,, is taken out of the circuit and the standard is inserted. The 
circuit is again balanced using arm C only, giving a resistance reading R, and 
a capacitance reading C,. Knowing the resistance of the coil under test and 
that of the standard we now have enough data to compute the effective induc- 
tance by use of formula (3). 

The resistors were tested at two frequencies, 60 and 500 cycles per second. 
At a frequency of 500 cycles per second measurements were made for a series 
of values of current through the shielding coil, on all three of the units. A 
number of measurements were made on the 10,000 ohm unit and also the 
5000 ohm unit at 60 cycles. Since the measurements at 60 cycles on the 5000 
ohm unit were obtained with difficulty due to low sensitivity, no attempt was 
made to obtain 60 cycle measurements on the 1000 ohm unit. 


8. RESULTS FOR IDEALLY SHIELDED RESISTORS 


The data obtained on the three units are listed in the tables below. Stand- 
ards with 20 cm spacing were used. The first column of the tables is headed 
I, the current through the working coil, and the second is headed by J,, the 
current through the shielding coil. The other columns are headed by symbols 
whose meaning has already been explained. 

The last column of the above three tables gives the effective inductance 
as observed for various currents through the shielding coil. These values are 
listed to the nearest hundredth millihenry, since computation to a greater 
degree of accuracy is not justified in measurements on resistances of the order 
of magnitude dealt with here. The individual figures set down in the table 


TABLE I. 10,000 ohm unit. 








Te l, R, C R 








w og C. Lw 

m.a. m.a. ohms mmf ohms mmf mh 

500 cycle measurements 
10.6 0 9930 77.9 9982 68.7 —1.44 
10.6 2.65 9930 80.2 9982 74.4 —1.10 
10.6 5.30 9930 74.8 9982 73.4 —0.67 
10.6 7.95 9930 68.8 9982 70.4 —0.37 
10.6 10.60 9930 69.2 9982 74.1 —0.04 
10.6 13.25 9931 57.0 9983 65.8 +0.3 
10.6 15.90 9931 52.7 9983 64.1 +0 .60 
10.6 18.55 9931 44.2 9983 60.8 +1.01 

60 cycle measurements 
10.6 7.95 9929 63.2 9982 63.8 —0.47 
10.6 10.60 9929 56.4 9982 61.0 —0.08 
10.6 13.25 9929 48.8 9983 58.7 +0.45 
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TABLE II. 5000 ohm unit. 



































Ie I, R, at x, Co | 
m.a. m.a ohms mmf ohms mmf mh 
500 cycle measurements 
12.0 0.0 9934 80.4 9970 76.8 —0.25 
12.0 3.0 9934 84.3 9970 81.2 —0.22 
12.0 6.0 9934 81.5 9970 80.2 —0.13 
12.0 9.0 9933 79.1 9970 79.3 —0.06 
12.0 12.0 9935 70.0 9973 71.4 0.00 
12.0 15.0 9936 65.3 9973 68.0 +0 .06 
12.0 18.0 9936 66.0 9973 70.3 +0.14 
12.0 20.0 9933 69.3 9970 a +0.22 
60 cycle measurements 
12.0 12.0 9932 60.5 9964 61.3 —0.04 
12.0 15. 9933 57.9 9964 59.7 +0.01 
TABLE III. 1000 ohm unit, 500 cycle measurements. 

Ie I, R, C. Re Co Le 
m.a. m.a. ohms mmf ohms mmf mh 
12.0 0.0 9981 64.5 9997 61.9 —0.01 
12.0 3.0 9981 64.5 9997 62.5 —0.01 
12.0 6.0 9981 62.5 9997 61.0 0.00 
12.0 9.0 9981 63.3 9997 62.1 0.00 
12.0 12.0 9981 61.5 9997 60.4 0.00 
12.0 15.0 9981 60.9 9997 59.7 0.00 
12.0 18.0 9981 63.3 9997 62.9 +0.01 
12.0 20.0 9981 63.0 9997 62.3 








represent the average of from five to ten separate measurements. The resist- 
ance balance changed only slightly, because of heating in the standard, during 
a series of measurements. The individual capacitance balances, however, 
were subject to some variations. In the case of the 10,000 ohm unit the aver- 
age deviation from the mean at 500 cycles for ten capacitance readings was 
approximately one percent. At 60 cycles, however, the average deviation 
from the mean for the same number of readings was seven percent. For the 
other two units this variation was still quite small at 500 cycles, not exceeding 
two percent; but at 60 cycles it mounted to about 10 percent for the 5000 
ohm unit. No 60 cycle measurements were attempted on the 1000 ohm unit, 
for the limit of sensitivity of the detector was approached in the measure- 
ments on the 5000 ohm unit. Since the 60 cycle measurements agreed closely 
with the 500 cycle measurements for the two higher resistance units, it is 
reasonable to assume that the same would have been true for the 1000 ohm 
unit. The sensitivity at 500 cycles was sufficient to measure time constants 
of the order of 10-8 seconds to a fair degree of accuracy, even for the 1000 
ohm unit where the results were consistent to within two percent. 

The curves of Fig. 4 illustrate the nature of the variation of the effective 
inductance with the variation of the current J,, through the shielding coil, 
for the three units. It is seen to be linear with J, as may be expected, and 
changes from negative to positive as J, is increased from zero. The slope of 
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the curves have been found to be directly proportional to the product of the 
distributed capacitance between the two coils and the resistance. Silsbee® 
found a similar result for resistors of the type illustrated by Fig. 1B. The 
point on each of the curves of Fig. 4 where the two currents are equal are 
marked by leaders. This point is the most important, even though the phase 
angle is not quite zero, because it represents the conditions under which the 
resistor is used in practice. 

The change of sign of the effective inductance is due to the reversal of 
the direction of capacitance currents between the two coils. When the current 
in the shielding coil is less than that in the working coil, the potential of any 
point on the shielding coil is lower than that of a corresponding point on the 


EFFECTIVE INDUCTANCE - MILLIHENRIES 





Fig. 4. Curves showing the variation of effective inductance with the variation of the current 
through the shielding coil for the three ideally shielded resistors. 


working coil. Therefore, the capacitance current passes from the working 
coil to the shielding coil. Now, when J, is greater than J,,, the potential of 
any point on R, is higher than that of a corresponding point on R,,, and thus 
the capacitance current passes from R, to R,. This reversal of capacitance 
currents corresponds to a phase change of 180° and hence accounts for the 
change of sign of the effective inductance. 

In Fig. 4 it will be noted that the effective inductance at the point J,=J, 
becomes increasingly more positive as the resistance of the unit decreases. 
It is still slightly negative for the 10,000 ohm unit but becomes decidedly 
positive for the 1000 ohm unit. It will be further noted that the agreement 
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between the results at the two frequencies, though close in all cases, becomes 
less close as the resistance of the unit decreases. This is explained by the 
fact that in the lower resistance units the effect of self-inductance becomes 
appreciable compared to the capacitance effects. In the 10,000 unit the self- 
inductance is extremely small compared to the capacitance and therefore the 
curves at the two frequencies almost coincide. 

Since the values of the phase angle with equal currents through the two 
coils are the ones in which we are most interested, they are set down in Table 
IV as computed for a frequency of 60 cycles. The agreement between the re- 
sults at the two frequencies was close enough to justify taking the simple 
average of the two sets of data where measurements were made at two fre- 
quencies. 




















TABLE IV. 
Resistor je Phase angle = 0 Tan 0=0 
10 ,000 ohm —0.06 mh —0.47 sec. —0.0000022 rad. 
5000 ohm —0.02 mh —0.31 sec. —0.0000015 rad. 


1000 ohm +0.005 mh +0.19 sec. +0 .0000019 rad. 











Referring to Table III it is seen that the effective inductance for equal 
currents in the two coils of 1000 ohm unit is set down as zero. This simply 
means that the value is less than 0.005 mh, since computations were made to 
the nearest 0.01 mh. Although the error in measurement was probably only 
a fraction of 0.01 mh, computation to a higher degree of accuracy was not 
necessary in this work. The purpose of the measurements on these shielded 
resistors was to show that the phase angle is less than 10~° radians, or, that 
the time constant is less than 10-7 seconds. 0.01 mh, even for the 1000 ohm 
unit, which is the worst case, corresponds to a time constant of 10~* seconds, 
which is negligible for all but the most precise measurements. By drawing an 
average straight line for the points obtained for the 1000 ohm unit, Fig. 4 
indicates that the effective inductance is decidedly positive and less than 
+0.005 mh for the point 7,=J/,. Thus, the value +0.005 mh set down in 
Table IV for the 1000 ohm resistor is probably larger than the actual value. 
This corresponds to a time constant of 5X10~-® seconds. 


9. EXTENSION OF THE RANGE FOR THIS METHOD 


In the preceding pages it has been shown quite conclusively that the cir- 
cuit gives consistent results for resistances ranging from 1000 ohms to 10,000 
ohms and can be used to measure time constants as small as 10~* seconds 
with a fair degree of accuracy. The type of shielded resistances upon which 
these measurements were made may be said to possess the lowest phase angle 
of any shielded resistor found in practice and thus the most difficult to meas- 
ure. Even more satisfactory results may be expected in measurements on 
other types of shielded resistors. By obvious and very simple modifications 
the circuit may be adapted to either of the other two types of shielded resis- 
tors illustrated in Fig. 1. 
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The upper limit of this circuit is not 10,000 ohms but has been used by 
the writer for resistances up to 50,000 ohms. Since parallel wire standards 
having a resistance of more than 10,000 ohms become quite long, measure- 
ments were carried out on the higher resistances in a slightly different man- 
ner, in which a number of 10,000 ohm shielded resistors were necessary. The 
effective inductances of severa! 10,000 ohm resistors are measured separately 
using the 10,000 ohm standard. Once the effective inductance is known the 
appropriate number of 10,000 ohm resistors can be connected in series and the 
combination used as a secondary standard replacing the primary parallel wire 
standard. 

Fig. 5 gives the connections for measurements on a 20,000 ohm shielded 
resistor of the type illustrated in Fig. 1B. Arms C, D, and W, are always left 


WV 























Fig. 5. Circuit for comparing the phase angles of resistors. 


the same as in Fig. 3. The resistance of arm B is made approximately equal 
to that of the unknown in arm A. W, is made one-tenth that of arm A. The 
two resistances in the upper part of arm A are previously measured with the 
10,000 ohm standard. The two are connected in series and the resultant ef- 
fective inductance of the series connection is the sum of those for the in- 
dividual units. With the series connection of units as a secondary standard, 
the effective inductance of the 20,000 ohm unit is measured in exactly the 
same manner as that for a 10,000 ohm unit. The shields of the various units 
in arm A are maintained at the proper potential by connecting them to W; 
as indicated in Fig. 5. It has been found that with this circuit it is immaterial 


where the ground is placed, but to be consistent it was always connected to 
the midpoint of W,. 
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By the use of a series connection of units the circuit has been used with 
good sensitivity both at 500 and 60 cycles for resistances as high as 50,000 
ohms. As the magnitude of the resistance increases, however, the voltage 
must be raised to get the required sensitivity. For 50,000 ohm resistors about 
500 volts had to be used in order to get good results at 60 cycles. 

The circuit can also be used to measure the phase angle of ordinary 
shielded resistance boxes. For precise measurements such as the measurement 
of power factor of dielectrics, the phase angle of the average resistance box 
is not always negligible. The effective inductance of a 10,000 ohm box made 
by a well-known manufacturer has been found to be — 2.54 millihenries. This 
corresponds to a phase angle of —0.33 min. at 60 cycles. Before making the 
measurement the box was shielded in a manner illustrated in Fig. 1A. 


10. SUMMARY 


The small, but not always negligible phase angles of “non-inductive” 
shielded resistors are often ignored because of the supposed extreme difficulty 
of their measurement. The substitution method described in detail can be 
carried out without the use of any very special apparatus with the exception 
that it must be possible to vary the resistance of the box in arm C of Fig. 3 
by as little as 0.01 in 10,000 ohms. The use of only one arm of the bridge 
network in balancing, and the use of an air capacitor to balance for phase, 
makes it a simple matter to eliminate the effect of parasitic influences. The air 
capacitor as used here is subject to no parasitic effects. The measurements on 
the ideally shielded resistors prove conclusively that consistent results can 
be obtained with an accuracy sufficient for the most precise measurements in 
which such resistors may be used. Time constants as low as 10~* seconds 
were measured with a fair degree of accuracy. 

The circuit is easily adaptable to all types of resistors. By the use of sev- 
eral resistors as secondary standards the range of the circuit can be extended 
to 50,000 ohms. The circuit can also be used to measure the phase angles of 
resistance boxes. The average resistance box has an inherent phase angle of 
sufficient magnitude to be of importance in many precise measurements. 

It may be said, therefore, that the circuit described herein affords a com- 
paratively simple means of measuring the phase angles of resistances ranging 
from 1000 ohms to 50,000 ohms. For frequencies as low as 60 cycles this range 
is narrowed slightly; it being from 5000 to 50,000 ohms. 
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NE OF the achievements of ‘modern physical science is the success with 
which the structure of atoms and molecules has been elucidated. It has 
been possible to assign definite structures to a number of the simpler mole- 
cules and to determine the molecular weights of an enormous number of more 
complicated ones. However, serious difficulties arise when it becomes neces- 
sary to associate a definite weight with the complex molecules or other struc- 
tural units which form such important substances as rubber, cellulose, body 
tissue, protoplasm, gelatin, glue, haemoglobin, egg albumin, synthetic resins, 
wool, silk, and the cellulose esters. The scientist and technician must now 
use these materials without having much information about their structures. 
In spite of these difficulties modern physical studies are revealing the size 
and weight of the more elementary units of the structures, and the manner 
in which they are arranged to form the massive substance. It is with the first 
of these problems that we are interested. A successful method to determine 
the micellar or particle weight of such substances consisting of spherical mole- 
cules or particles has been developed by Svedberg and his associates. It con- 
sists of a study of the velocity with which a particle whose weight or size is 
to be determined sediments in a centrifugal field of known force. The tend- 
ency of the particle to settle out in such a field is opposed by the linear fric- 
tional resistance offered to its movement in a viscous medium in which it is 
suspended. Well-known laws enable us to calculate the size of a sedimenting 
particle since the force of the centrifuge just overcomes the force of friction 
when the particle falls with a constant and measurable velocity. 

Another method by which the size and therefore the weight of a spherical 
particle can be determined requires experimental data of quite another type 
but at the same time embodies certain similar principles. It has recently been 
described by Marinesco,' and independently but somewhat later by Williams 
and Oncley.?” It is an application of the dipole theory of Debye and depends 
upon the fact that because of the frictional resistance of the medium to the 
rotation of a suspended particle a finite time is required for its orientation in 
an alternating electrical field. According to this theory it follows that if the 
frequency dependence of the dielectric constant is determined for a system 
composed of electrically dissymmetrical particles suspended in a liquid me- 


1 (a) Marinesco, Comptes Rendus 189, 1274 (1929); see also Conference-communication 
faite a la Soc. de Chimie-physique de France. Seance du 26 juin 1929; (b) J. chim. Phys. 28, 
51 (1931); (c) Kolloid Zeits. 58, 285 (1932). 

2 (a) Williams, J. Frank. Inst. 211, 581 (1931); (b) Williams and Oncley, J. Rheology 3, 
271 (1931). 
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dium there will be found to be a region in which the dielectric constant de- 
creases as the frequency is increased. The size and therefore the weight of the 
suspended particle is expressed in terms of the dielectric constants for zero 
and for infinite frequencies, and either the time of relaxation or the true inner 
frictional coefficient of the liquid. 

The qualitative argument may be briefly stated by saying that at suffi- 
ciently high frequencies the dielectric constant will fall off, because the inner 
friction constant of the medium prevents the orientation of the dipole par- 
ticles which contribute to the dielectric constant. The polarization consists 
of two parts, a polarization due to deformation of the particle by an adjust- 
ment of the atoms and electrons, and a polarization due to the orientation of 
the particle itself. The transition from both distortion and orientation effects 
to the distortion effect alone will occur in a frequency region defined ap- 
proximately by the equation 


yr = | (1) 


where 7, called the time of relaxation of the molecules, is the time required 
for 1/e of the particles to assume a random distribution after the applied field 
has been removed, and »y, is the frequency of the alternating field. It should 
be emphasized that this equation serves only to indicate the region in which 
anomalous dispersion occurs, and means simply that v, and 1/7 must be of 
the same order of magnitude. 

In order to have the situation clearly before us it seems worth while to 
give an abbreviated outline of the exact dispersion theory as it has been given 
by Debye, and then to extend the theory so that it may be applied to solu- 
tions or suspensions of polar particles in dilute solution in a nonpolar solvent. 
It will be seen that as a result of these calculations the transition between the 
high and low values of the dielectric constant will take place in the frequency 
region defined by the expression, 


0.1 


IA 


x <= 10. 
In this equation 


€9 +2 xrtv 
€é. + 2 kT 


where €, is the dielectric constant at zero frequency, e€,, is the dielectric con- 
stant at infinite frequency, k is the Boltzmann constant, 7° is the absolute 
temperature, v is the electrical frequency in cycles per second, ¢ is the “inner- 
frictional constant” which, when and if Stokes’ law can be applied is equal to 
82rr*n, n is the ordinary coefficient of viscosity of the solution and r is the 
radius of the particle. 

If the system is a polydisperse one the curves obtained when either of the 
dielectric constants (real or imaginary) or the power factor are plotted against 
the frequency will be flattened out. They may be regarded as the result ob- 
tained when the simple curves for the several particle sizes are added to- 
gether. 
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The treatment of Debye’ resolves itself into a consideration of the gen- 
eralized distribution of the dipole moments of the polar molecules if affected 
by an alternating electrical field. Saturation effects do not enter in because 
the field strengths used in the ordinary experiments are small. The case of a 
field of force in one direction only in which the field is variable with time but 
constant in direction is considered sufficient for the dispersion problem. If 
we now suppose a spherical surface drawn in the liquid and take the center 
of this sphere as the pole of our polar coordinate system, then obviously the 
distribution function for the dipole moments must depend upon the polar 
angle @ and upon the time ¢. Thus the number of molecules whose moments 
point in the direction of a solid angle dQ at a given time can be denoted by 
fdQ where f is f(9, ¢). 

Let us consider the events taking place during an interval of time 6¢. If 
t=0 at the beginning of this interval at the end of it we will have ¢=6¢. Let 
us further assume that 6¢ is sufficiently large for all molecules whose moments 
were in dQ. at t=0 to have shifted outside this angle when ¢ = 6/, and yet small 
enough so that during the time 6¢ the moments will shift not more than a few 
degrees. Then the number of molecules whose moments enter dQ in the time 
6t will be 

5t(0f/dt)dQ = Ay + Ae (2) 


where A, is the contribution due to rotations produced by the impressed field, 
and A, is the contribution due to rotations caused by the Brownian move- 
ment. 


The result of the calculations for A; and A, will be given at once: 











A : = ( . it si ) dQ (3) 
=- —{ f-—-ét sin 
, sin @ 006 f ¢ 
©? /cos@ a 0? 
a, = aa—(“ af *, (4) 
4 \sin 6 00 op? 


The negative sign in the expression for A, takes account of the fact that the 
torque acts opposite to the direction of increasing angle @. 

In these equations the torque on a molecule is JJ = —yF sin 6, where yp is 
the electric moment and F is the electrical intensity acting upon one molecule 
or particle. The significance of © will be understood from the fact that the 
constant ©?/6¢ is closely related to the well-known constant £2/6¢ in the ordi- 
nary expression for the Brownian movement. It is the mean angular rotation 
squared per unit interval of time. 

Substituting Eqs. (3) and (4) in (2) there results the expression 


af . ae ©? af M 
St a e 
dt sin @ 30 4at 00 ¢ 


The constant 0?/45¢ evaluates to kT, ¢, thus, 





3 Debye, Polar Molecules, Chemical Catalog Co., New York, 1929. German Edition, Hirzel, 
Leipzig, 1929. 
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0 1 0 0 
¢ Z oa sin o( er J _ u-f)| (6) 
ot sin 6 06 00 ; 


which is the differential equation expressing f as a function of @ and ¢. 

The intensity of an electrical field which is alternating periodically is 
given by the formula F= Fp) cos wt where w is the number of vibrations in 
2m seconds (vibrations in one second). For mathematical simplicity this is 


expressed as the real part of F= Fye”’'. Therefore, \J = —pFye™! sin 0. Making 
the further substitution r={/2k7 Eq. (6) becomes 
of 1 0 of ful’ oe'*' 
2— = ——- sin of —~+$- = sin 0)|. (6a) 
dt = sin @ 06 00 kT 
The solution of this equation is 
1 pl oe’! 
f=A E + —.- —— COS o|. (7) 
1+ wr kT 
The mean moment of a particle or molecule is 
ph? 1 
m = (as + meter —)F. (8) 
3RT 1 + twr 


The meaning of the complex moment is that there is a difference in phase 
between the dipole moment and the internal force. 

Also, since m=aF and (€—1)/(€+2) = (4r/3)na where a is the polariza- 
bility and is the number of polar particles per cubic centimeter, we may 
write 


e—1 4r pe? 1 d 
cuneate he | as +—-—- —| = Pw) —- (9) 
e+2 3 3kT 1+ iwr G 
Thus, the molar polarization P(w) and dielectric constant ¢ for a polar 
substance of deformation polarizability ao, electric moment yp, density d, and 
molecular or particle weight G has been calculated assuming that the Clau- 
sius-Mosotti hypothesis is valid for the pure substance. However, this is 
known to be the case only when such polar molecules are dissolved or other- 
wise suspended in a nonpolar solvent medium. Accordingly it is necessary 
to give consideration to the adaption of these equations to such a system. 
The following is such ‘a treatment for a solution in which m2 polar molecules 
per cubic centimeter whose deformation polarization is a") and whose dipole 
moment is pe are dissolved in m, molecules per cubic centimeter whose def- 
ormation polarization is a;“, and whose dipole moment is zero. 
Thus, Eq. (9) above becomes 


-—— 1 4n Nope” 1 
———_ = — |] nayio)t Nea) + --——-— | (9a) 
€i2 + 2 34 3kRT 1+ twre 
d 
= Pi(w) —— 
la * fle 
€2(2%) — 1 1 (a — 1 €2(%) — ) 
7 €:2.(%) + 2 1 + iwr\e2(0) +2 en(~) +2 
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Therefore 
’ ae 1 - ba (fiGi + fG2) 
q13 = 19 — 33 = 


€12(% — 1 1 €12 0) —1 €12 x) — 1 
1+2[! ) p(s = — ax(*) boa 


ex(%) + 2 1 + wre €12(0) re 2 €12(2) +2 


























—_— 4 a — (10) 
OR POP) 
€2(©) +2 1 +iwre \e2(0) + 2 €9(%) + 2/] 
7 €12(0) + xei2(2 ) 7 £1200) — eas) (10) 
1 + x? 1+ x? 
where 
€:2(0) + 2 
a 
€12(% ) oh 7 
Thus, 
(0 = 9( © 
€12) = €12(% 2) +— is (11) 
and 
0) — 
‘te at 


1+ x? 


In these equations the dielectric constant for the very low frequencies is 
designated as €9 for the pure substance and as «2(0) for the binary mixture. 
The quantities ¢,, and €2(%) represent the corresponding dielectric constants 
for very high frequencies. The molar polarization is designated as P(w) in 
order to distinguish it from that of the static case. The quantities f; and fe 
are the mol fractions of the two components of the system. 

With these equations there may be calculated the expected behavior of 
any binary system in which the assumptions underlying both the law of 
Stokes and the Clausius-Mosotti hypothesis are fulfilled. It is necessary to 
obtain dielectric constant and density data for dilute solutions of the polar 
molecule or particle in a nonpolar solvent, and then extrapolate the results 
to zero concentration in order to find the behavior of the solute molecule un- 
der the assumed conditions. In making this extrapolation it is difficult to say 
just which value for the coefficient of viscosity 7 should be used. This extrapo- 
lation leads to additional difficulties if the law of Stokes cannot be applied. 
A detailed description of these difficulties involved in the use of the ordinary 
coefficient of viscosity and Stokes’ law for calculating molecular sizes from 
dispersion data has already been given by us”® so it is not necessary to discuss 
them at this point. 

The first application of the dipole theory to the study of colloidal chem- 
istry was probably made by Errera.** In 1928 Wo. Ostwald*® began the in- 
vestigation of the relationship of dielectric constant, polarization and electric 


‘@ Errera, Kolloid Zeits. 32, 157, 240, 273 (1923). 
* Wo. Ostwald, Kolloid Zeits. 45, 56 (1928), and a number of later interesting articles. 
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moment to the stability of a large number of dispersed systems. Shortly after- 
ward Marinesco' recognized that the dipole theory could be used to de- 
termine micellar size, (grandeur), and to explain a number of peculiarities 
observed by him during the course of his previous studies of the dielectric 
constant of various gelatin solutions. In the short note published in Comptes 
Rendus he states that “the anomalous dispersion of the dielectric constant 
of a colloidal solution whose micells have permanent moments permits the 
determination of the molecular volume, inasmuch as the time of relaxation 
of a micell is given by 


Tt = Aanor®/ kT => Xo/ r 


where 7 is the coefficient of viscosity, ¢ is the velocity of light and Xo is the 
wave-length corresponding to the beginning of the fall of the dielectric constant 
of the solution toward the region of the anomalous dispersion, and is deter- 
mined by the measuring circuit. The molecular volume in solution is 


where R is the gas constant and M is the molecular weight.” By defining Xo 
as the wave-length corresponding to the beginning of the fall of the dielectric 
constant of the solution Marinesco has clearly recognized the fact that a 
quantitative result cannot be obtained using the expression v.r=1 or any 
equation equivalent to it. The point at which the dielectric constant begins 
to fall will certainly be difficult to locate with any degree of exactness. It may 
best be done by the use of the formula for the dielectric constant which marks 
the beginning of dispersion, 


€o + 4re,,[(eo + 2)/(e, + 2)]? 


1 + 4n?[(eo + 2)/(e, + 2)]? 





In his more recent publications!» Marinesco has apparently considered 
it unnecessary to apply the exact relationships. Thus in the first of these ar- 
ticles!» the expression for the molecular weight of a colloid is given as, 


and Xx is clearly indicated in an accompanying figure to be the wave-length 
at which 


This is equivalent to the assumption that v.r=1, discussed above. In the de- 
termination of the molecular weight of gelatin Marinesco has worked with a 
0.1 percent solution at py=3.5. At \=6.5m, the dielectric constant was 
found to be 86.0, while at \ = 1.9m the value was less than that of pure water. 
Experiments at \=6.2m showed that the dielectric constant value had de- 
creased to 82.5. Assuming with him that the “region of anomalous dispersion 
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is consequently ‘very nearly’ the wave-length \=6m” the molecular weight 
is found to be 11,300. A previous determination of v, a quantity necessary for 
the calculation, had given 1.294 cm*. This result is compared with one re- 
ported by Cohn,’ M=10,300, which had been obtained in solutions of like 
concentration, but using other methods. 

The second of these articles, whose purpose is to review previous work, 
contains equally definite statements with respect to which objections may be 
raised. Thus it is stated that the dielectric constant curve passes through a 
point of inflection (Wendepunkt) for a wave-length \¢ corresponding to the 
time of relaxation 7, of the dissolved molecule, and that in this instant the 
period of the applied field is equal to 7,; thus 

ho /E = Arr no/ kT 
and 


1 do’ RT 
a € UNo 


M 


II 


In spite of the fact Marinesco did not trace the exact course of the dispersion 
of the dielectric constant, Ao’ is now assigned the value 6m, and the gelatine 
is again found to have the molecular weight 11,300. 

But it can be seen from Eqs. (11) and (12) that the real dielectric constant 
e’ will take on a mean value and the imaginary dielectric constant e’’ will 
pass through a maximum at the point where x =1. Under these conditions 


é9+2 xv, 
«+2 kT 
In other words the critical frequency y, is given by the relation 
€.+2 kT e +2 1 
eo +2 mf  e+2 2ar 
€. + 2 





= 1. (13) 
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Eq. (13) is the exact equation which should always be used in work of this 
sort. The equation for the molecular or particle weight should therefore be 
written in the form 
é~+2 rX RT 

—-—-- (14) 
é9 +2 2ré nv 








1 
M =— 
3 


If one recalculates the result of Marinesco for the particle weight of gela- 
tin (dielectric constant data from dilute solution at py =3.5) using our for- 
mula (14) it is found that the value may be of the order of magnitude 1/10 
of that previously reported.' Thus the excellent agreement indicated between 
the result obtained from the dipole theory and those obtained in other ways*® 
is apparent rather than real. 


6 Cohn, J. Biol. Chem. 63, 721 (1925). 
* Cohn, J. Biol. Chem. 63, 721 (1925); Biltz, Zeits. f. phys. Chem. 91, 705 (1916). 
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This difficulty of a low molecular or particle weight resulting from the 
application of the dipole theory is to be expected if results published from 
other laboratories are considered. To mention a single example Wyman’ has 
made much more complete and more exact dielectric constant measurements 
on solutions of the protein zein dissolved in 70 percent n-propyl alcohol over 
a temperature interval from 0° to 80° and at wave-lengths from 4 to 260 
meters. The data show the presence of anomalous dispersion in a convincing 
manner. In discussing the question of the relaxation time Wyman states “It 
was calculated on the basis of Stokes’ formula that the relaxation time of a 
spherical molecule of the size of zein (gram molecular weight 100,000) in a 
medium of viscosity 0.01 ought to be 8X10~* sec. The relaxation times, 
2.8 10-* and 5.4 X 10-8 just computed by the Debye theory from the experi- 
mental results both correspond to viscosities considerably greater than this, 
but since relaxation time and viscosity are proportional to one another they 
may be reduced for comparison. By multiplying 2.8 X 10~* by 0.010/0.017 and 
5.4X10-8 by 0.010/0.034 (the fractions being appropriate viscosity ratios) 
we obtain respectively 1.7X10-* and 1.610~*.” Both of these values are 
about { of that calculated by the theory by assuming Gz= 100,000 and using 
the proper mathematical relationship, which means that any molecular 
weight calculated from these data will be too small by the same factor. The 
data themselves seem to be at least as good as any which have been reported 
to date for a protein in a solvent. 

Two outstanding difficulties which must be overcome before Eq. (14) may 
be applied to the determination of a particle size or weight are given below. 
Because of them and because of lack of knowledge of the proper shape factor 
to use in connection with the orienting particle a really quantitative result 
has not as yet been obtained for any system in which the frequency variation 
of the dielectric constant is involved. When and only when these difficulties 
can be overcome will it be possible to determine a particle size or weight of 
true significance. 

1. The ordinary macroscopic coefficient of viscosity does not seem to 
measure the resistance to the orientation of the suspended or dissolved par- 
ticles. Instead there must be used what might be termed a microscopic vis- 
cosity. This is especially true when substances having larger molecules are 
used as solvent media for smaller molecules. 

2. The development of the theory is such that nonpolar suspension or 
solution media are required for a quantitative interpretation of the results 
either in terms of the size or electric moment of the suspended particle. Thus 
the application of the Debye theory to experiments made in polar solvent 
media will always be of highly questionable significance. 

It is interesting in this connection that Wyman extended his measure- 
ments to the gels which formed when the zein systems were allowed to stand. 
No change was observed in the dielectric constant as the gelation took place. 
The suggestion was made by him that the effective viscosity of the gels with 


7 Wyman, J. Biol. Chem. 90, 443 (1931). 
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respect to the rotation of the particles is not the ordinary viscosity entering 
the simple Stokes’ expression. This failure of the dielectric constant to change 
on gelation reminds one of the corresponding failure of the electrical conduct- 
ance of a salt solution containing gelatin to change to an appreciable extent 
when the system sets. 

The difficulty involved in the determination of the proper viscosity to use 
in the formula for the time of relaxation, and therefore for the molecular 
weight, has already been mentioned. We have even suggested that the theory 
might well be used to determine the true inner friction constant of a liquid, 
providing care is taken in the choice of system used. The system studied must 
conform to the assumptions which have been made in the development of the 
Clausius-Mosotti law, that is, the second difficulty mentioned above. 

Recent experiments with systems in which both rosin oil and castor oil 
molecules have been dissolved in various hydrocarbons indicate that the time 
of relaxation for these molecules may be practically independent of the mac- 
roscopic viscosity of the solution or solvent as measured with a rotating cyl- 
inder viscometer. In other cases where vaseline-like substances have been 
added to increase the viscosity of the solution anomalous effects have been 
observed. These results will be presented in another place. 

The frequency-dielectric constant curves make it possible to distinguish 
quite readily between a monodisperse and a polydisperse system. In the case 
of the polydisperse system the dielectric constant change is a very gradual 
one, spread out over a considerable range of frequencies, but for the mono- 
disperse system the change from the high value to the low value should be 
quite sharp. Recent work of Kitchin* on rubber systems gives definite indica- 
tion of an electrical inhomogeneity, which may or may not be accounted for 
on the basis of a dipole mechanism. The same thing is true of work by Race® 
on the high-frequency electrical properties of insulating oils, in which the 
change in the dielectric constant appears to represent the sum of the contribu- 
tions of particles of several sizes. The average size of the particle radius as 
calculated by Race is 3.4X10-* cm, a value which must be too small by a 
considerable factor. It may be remarked that a number of the suspended 
particles in these oils can be seen to scatter light when observed in the ultra- 
microscope. Such particles should have radii at least eight to ten times as 
large. It may also be mentioned that no correlation was observed between the 
change in the viscosity of the oils and the frequency region in which the 
anomalous dispersion took place. 


SUMMARY 


The dipole theory of Debye has made it possible not only to learn some- 
thing about the structures of simpler molecules but also to approximate the 
size and therefore the weight of a dissolved or suspended molecule or particle. 
The latter result should follow, as has been pointed out by Marinesco, and 


§ Kitchin, J. Ind. Eng. Chem. 24, 549 (1932). 
* Race, Phys. Rev. 37, 430 (1931); J. Phys. Chem. 36, 1928 (1932). 
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also by us, from a study of the time of molecular rotation in an alternating 
electric field. The method depends upon the fact that because of the frictional 
resistance of the medium to the rotation of the suspended particles a finite 
time is required for their orientation in the field. If the frequency dependence 
of the dielectric constant is determined for a system composed of electrically 
dissymmetrical particles suspended in a liquid medium there will be found to 
be a region in which the dielectric constant decreases as the frequency is in- 
creased. The theory expresses the size and therefore the weight of the sus- 
pended particle in terms of the observed critical frequency and the “viscos- 
ity” of the medium. 

It is shown that the theory as used by Marinesco can give only an approxi- 
mate result for a molecular weight. The exact theory of Debye is developed 
only to the extent which we believe to be necessary to demonstrate this fact. 
In addition we have considered the changes necessary to the treatment of 
Debye to make it applicable to the type of binary mixture which must be 
used in the experimental work. The data now available in the literature are 
as yet insufficient for the exact determination of the particle size because the 
systems involved fail to conform to the assumptions underlying either the 
law of Stokes or the Clausius-Mosotti relation. Nevertheless, calculations 
from them indicate that once it becomes possible to evaluate the true inner 
friction constant of the system important results may be obtained. The 
method readily distinguishes monodisperse and polydisperse sols. 


Note added in proof: Since this article was submitted, Bloch and Errera 
have published an article in Phys. Zeits. 33, 757 (1932) in which they come 
to much the same conclusions that we have. 
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Capillary Action in Impregnated Paper Insulation 
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The capillary rise of oil in vertical strips of impregnated paper obeys the same law 
as the rise of liquids in capillary tubes of circular cross section. From the viscosity 
and surface tension of an oil and its rate of rise in a given paper, it is possible to com- 
pute the “effective capillary radius” of the capillary pores of the paper. The rate of 
rise in a calibrated paper affords a convenient means of measuring the surface tension 
when the viscosity of the liquid is known. The rate of penetration of an oil into a paper 
is directly proportional to the “penetrativity” of the oil as defined in the theory of 
capillary action. The penetrativities of nine insulating oils and the effective capillary 
radii of seven different papers have been determined. A uniform relationship be- 
tween the effective capillary radius and the Gurley air resistance of the paper is in- 
dicated. 


N THE course of a study of the electrical properties of impregnated paper 

as related to the impregnating process, a series of observations was made 
of the rate of capillary rise of a number of oils in vertical strips of several 
types of paper. It was found in every case that the rise was accurately pro- 
portional to the square root of the time, over the interval studied, 1 hour to 
100 hours. 

The paper strips, 2.54 cm wide, were mounted vertically, with their lower 
ends dipping in the oils, and placed in a constant temperature box provided 
with a glass window. The progressive rise of the level of oil saturation was 
read with a cathetometer. Fig. 1 shows the results when using a well-known 
wood pulp paper, 0.004” thick, as used for cable insulation, and a number of 
different oils. Some of the constants of the oils are given in Table I. Fig. 2 
shows the results on a number of different papers (see Table II) with one 
type of oil. 

TABLE I. Physical properties of cable oils. August, 1932. 











Penetrative 
a Conduc- 
Pour Flash | Specific | Viscosity in poises ees ot = tivity 
Oil point point | gravity 7 — Fait Paper | mhos/cc 
°C °C 40°C 40°C 100°C dynes/cm cm-sec. 10-4 
40°C x 10-3 40°C 
40°C 
100 35.0 274 0.8788 7.5 0.235 51.8 S.34 34.2 
101 15.7 296 | 0.881 5.7 0.258 31.2 5.0 218.0 
102 —7.0 288 | 0.9268 7.2 0.27 20.1 3.55 1065 .0 
104 0.0 0.83 0.132 | 0.0295 | 32.0 33.0 1.0 
105 —7.0 282 | 0.8829 | 4.9 0.265 31.3 5.47 84.8 
106 4.4 0.935 417.5 0.446 21.6 2.28 149.0 
107 4.4 263 | 0.903 4.07 0.22 38.3 6.52 8.08 
108 —29.0 148 | 0.8805 | 0.1602 | 0.027 28.2 28.2 24.8 
109 —6.7 271 | 0.9015 2.46 0.155 33.0 7.67 22.0 
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TABLE II. Physical properties of cable papers. June, 1932. 








Air 
resistance App. 
Kind of paper Gurley specific 
seconds gravity 
8 mil low density 100 0.82 
8 mil high density 180 0.88 
5 mil high density 290 0.90 
4 mil kraft 640 
8 mil super calendered 640 1.07 
6.5 mil super calendered 1300 1.13 
5 mil super calendered 1.13 





Effective Penetrative 
capillary power 
radius of oil No. 109 
cm-sec. X 107% cm-sec. X 107% 
at 26.4°C 26.4°C. 
1.072 7.63 
0.987 6.25 
0.9325 5.77 
0.816 4.42 
0.783 4.07 
0.7465 be 
0.6975 3.235 
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Fig. 1. Penetration study. Kraft paper immersed in various oils. Height-time curves at 40°C, 
March, 1932. /=Kt'/?. K =penetrative power = (r/2)'/2- (y/n)"/?. 


The uniform relation: 
l= 


Ki'/2 








(1) 


shown in these curves is the same as that found for the rise of liquids in small 
capillary tubes of circular cross section. This relation is also accounted for in 
the theory of capillary action. For a straight vertical tube of radius r, Wash- 


burn! gives 


1 Washburn, Phys. Rev. 17, 273 (1921). 
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1? = [(y/n)(cos 6/2) |rt (2) 
(dl/dt) = (1/n)(y/4l) cos 8 (3) 


in which /=height of rise, y=surface tension, 7 = viscosity of the liquid, and 
@ is the angle between the inner surface of the tube and the meniscus of the 
liquid (6 =0 if the liquid wets the tube). It is also shown by Washburn that 
for very small values of r, 8 may always be taken as 0. In these expressions 
the negative (downward) hydrostatic pressure due to the small values of / is 
neglected as compared with the capillary (upward) force (=(2y/r) cos 6). 


and 
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Fig. 2. Penetration study. Various papers immersed in oil. No. 109. Height-time curves at 
26.4°C, June, 1932. /=K#/?, K =penetrative power=(r/2)"/2- (y/n)"/2. 


The quantity [(y/m) (cos @)/2) | is called the “penetrativity” of the liquid; 
it is seen to be equal to the square of the distance the liquid rises in a tube 
of unit radius in unit time. For convenience we will use the constant K of 
formula (1) in the discussion of our results, and will call it the “penetrating 
power.” The penetrating power involves the radius of the tube as well as the 
properties of the liquid, and so refers to the penetration of a particular liquid 
into a particular tube or paper. 

We note then that the expression / = K?!/* applies to a single fine straight 
capillary tube and also to the mesh of fine capillaries making up the paper 
strip. For the tube, the penetrating power, 


K = [(yr/n)(cos 0/2) ]*/? (4) 


or, for small tubes, 6=0, 
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yr\ 1/2 r\l/2 y\1/2 
, «fs =| — . pe A (5) 
ae “ayo 


In these values of K the only factor peculiar to the tube is its radius 7, y 
and 7 being the surface tension and the viscosity, respectively, of the im- 
pregnating liquid. 

Obviously, owing to the complexities of structure of the cellulose fibers in 
a paper tape, it is not permissible to assume that the various factors governing 
the rate of rise of a liquid in a vertical strip enter in the same way as in the 
case of a capillary tube. Nevertheless, the identity of the experimental rela- 
tion in the two cases is very striking, and it is immediately obvious that with 
a given liquid and paper there is a particular value of r for each paper, which 
enables us to compute the rise of an impregnating liquid in the paper by using 
the same Eq. (2), which applies to a straight capillary tube. We may call 
this value of r the “effective capillary radius” of the paper. The value of the 
effective capillary radius is readily computed from the values of K as found 
experimentally, if the values of y and 7 are known. 

The values of 9, the viscosity, are usually available or readily measured 
for any insulating oil. However the values of y, the surface tension, are by 
no means so easy to obtain, doubtless owing to the difficulty of accurate 
measurement for liquids of the high values of viscosity pertaining to insulat- 
ing oils: for viscosities above 2 poises, static methods are subject to serious 
error. We have met this difficulty by making direct measurements of y and 
n on two of the lightest oils studied and also on dibutyl phthalate. The last 
named liquid is practically non-volatile under ordinary conditions (hence is 
closely similar to the oils), and the values of y and 7 are given in the Interna- 
tional Critical Tables; our own measurements (Saybolt apparatus for 7, and 
Jolly balance for y) were in good agreement. In this way we were able to 
compute the effective capillary radius for any paper with three different liq- 
uids, and so “calibrate” the paper for the measurement of y for any oil for 
which 7 was known. In the case of a low viscosity insulating oil, the value of 
surface tension as measured by the penetration of the oil into a calibrated 
paper agreed within 2 percent with the value obtained by a series of measure- 
ments on a precision Jolly balance. 

We thus have a convenient indirect method for measuring the surface 
tension of the heavier viscous liquids, in the rate of rise in a capillary tube 
or in a calibrated paper. This method is so simple that it appears probable 
that others have used it. 

Table III gives the values of surface tension and of viscosity as measured 
at three temperatures for dibutyl phthalate, and for oils No. 104 and No. 
108, both thin water-white oils of the type used in oil-filled cables. The table 
also gives the values of K (cm and sec.) and the corresponding computed 
values for 7 the effective capillary radius of the 4 mil kraft paper. It will be 
noted that 7 is of the order of magnitude 10-5 cm, and that the value is fairly 
constant as based on measurements with different oils and different lengths 
of the same paper. However, r increases markedly with increase of tempera- 
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TABLE III. Capillary calibration of 4 mil kraft paper. June, 1932. 














Temperature of the Measurement 
Property Comp. No. 104 No. 108 Dibutyl phthalate 
24°C 40°C 60°C 24°C 23°C 40°C 60°C 

Specific gravity 0.84 0.83 0.818 0.891 1.044 1.0307 1.015 
Viscosity () in poises 0.26 0.132 0.073 | 0.34 0.178 0.0905 0.054 
Surface tension dynes/cm 

(y) 33.45 32.0 30.27 | 30.0 36.9 35.5 33.5 
Penetrative power cm-sec. 

x 10-*(K) 23.7 33.0 45.8 18.83 | 28.85 41.2 50.3 
Eff. capillary rad. cm 

X 10-5 (r) 0.875 0.899 1.011 | 0.804] 0.801 0.867 0.821 























TABLE III. (Continued). Effective capillary radius of 4 mil kraft paper.* 


24°C 40°C 60°C 
0.81610 cm 0.882 X10 cm 0.916 10% cm 


* Note: Computed from average values of the penetrative power of three different liquids 


into4 mil kraft paper. 
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Fig. 3. Effective capillary radius vs. temperature. Seven 
different papers, July, 1932. 


ture presumably due to volumetric expansion. Fig. 3 shows the effect of tem- 
perature upon the value of the effective capillary radius for several different 
papers. 

Table II gives the values of the effective capillary radius for a number of 
different grades of paper. The papers are arranged in the order of increasing 
air resistance. The influence of increased density and of super-calendering on 
the air resistance is readily seen. The fourth column of Table II, gives the 
effective capillary radius, computed as indicated above, from measurements 
of the capillary rise of oil No. 109. It is seen that the effective capillary radius 
decreases uniformly through the entire group, with increasing air resistance. 








as 














IMPREGNATED PAPER INSULATION 329 


The fifth column of the table gives the value of the penetrating power of oil 
No. 109 into the various papers and the same uniform correlation with the 
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Fig. 4. Effective capillary radius vs. Gurley air resistance 26.4°C, July, 1932. 
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air resistance is evident. The relation between air resistance and effective 
capillary radius is shown graphically in rectangular coordinates in Fig. 4, and 
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in Fig. 5 the same data are plotted in double logarithmic coordinates. The 
straight line of the double logarithmic record indicates the Gurley air re- 
sistance of the paper is inversely proportional to the 6.5th power of the ef- 
fective capillary radius. There is no obvious explanation of this relationship. 
If the longitudinal capillary channels were of the same type as the transverse 
air channels, the inverse air resistance should be proportional to the square 
of the effective capillary radius. The appearance of the 6.5th power makes it 
evident that the two types of channel are widely different. A wide variation 
in air resistance is accompanied by a much smaller variation in capillary 
radius. Increased density and super-calendering have a much greater influ- 
ence on air resistance than on capillary radius. Apparently the capillary 
properties arise in the narrow passages at the surfaces of contact of the cellu- 
lose fibers all compressed in the plane of the paper, while the air resistance 
obviously depends on the transverse density of the mesh as well as on the 
surface conditions. The two types of passage are obviously different, and 
were it not for the definite relation suggested by Fig. 5, one might readily 
conclude that the capillary and impregnating properties of a paper are more 
or less independent of the air resistance. That some connection exists how- 
ever, seems to be definitely indicated by the clean cut relation of Fig. 5. 

Further studies are under way as to possible relationships between the 
penetrating power of the oil, the capillary properties of the paper, and the 
electrical properties and behavior of the impregnated paper. 
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